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We give a surface integral derivation of post-l-Newtonian translational equations of motion for a 
system of arbitrarily structured bodies, including the coupling to all the bodies' mass and current 
multipole moments. The derivation requires only that the post-l-Newtonian vacuum field equations 
are satisfied in weak-field regions between the bodies; the bodies' internal gravity can be arbitrarily 
strong. In particular black holes are not excluded. The derivation extends previous results due to 
Damour, Soffel and Xu (DSX) for weakly self-gravitating bodies in which the post-l-Newtonian field 
equations are satisfied everywhere. The derivation consists of a number of steps: (i) The definition 
of each body's current and mass multipole moments and center-of-mass worldline in terms of the 
behavior of the metric in a weak-field region surrounding the body, (ii) The definition for each 
body of a set of gravitoelectric and gravitomagnetic tidal moments that act on that body, again in 
terms of the behavior of the metric in a weak-field region surrounding the body. For the special case 
of weakly self-gravitating bodies, our definitions of these multipole and tidal moments agree with 
definitions given previously by DSX. (iii) The derivation of a formula, for any given body, of the 
second time derivative of its mass dipole moment in terms of its other multipole and tidal moments 
and their time derivatives. This formula was obtained previously by DSX for weakly self-gravitating 
bodies, (iv) A derivation of the relation between the tidal moments acting on each body and the 
multipole moments and center-of-mass worldlines of all the other bodies. A formalism to compute 
this relation was developed by DSX; we simplify their formalism and compute the relation explicitly, 
(v) The deduction from the previous steps of the explicit translational equations of motion, whose 
form has not been previously derived. 



I. INTRODUCTION AND SUMMARY 

A. Background and motivation 

For slow motion sources in weak gravitational fields, 
general relativity can be accurately described in terms of 
a post-Newtonian approximation scheme. This approxi- 
mation scheme is extremely useful in applications and is 
very well developed. Reviews of post-Newtonian theory 
can be found in Refs. T, .2, 3] and in the textbook by 
Will 0. 

There are several different types of equations that arise 
in post-Newtonian theory. First, one has continuum field 
equations, which are usually specialized to gravity cou- 
pled to perfect or imperfect fluids. These have been 
derived up to post-2. 5-Newtonian order j^. At post-l- 
Newtonian order, they have been extended beyond gen- 
eral relativity to encompass the class of theories of gravity 
described by the parameterized post-Newtonian frame- 
work 4]. 

A second type of equation of motion applies to sys- 
tems consisting of N interacting, extended bodies mov- 
ing under their mutual gravitational interactions, in the 
limit where the bodies' sizes are small compared to their 
mutual separations. For such systems one has "point 
particle" equations of motion. Such equations were first 
derived Jl| at post-l-Newtonian order by Lorentz and 
Droste Q, and later independently by Einstein, Infeld 
and Hoffmann (EIH) 0. They were also independently 
derived by Petrova Q using a method devised by Fock 
0. These equations are usually called the EIH equations. 
In recent years the advent of gravitational wave astron- 



omy lid lll| has spurred renewed interest in such equa- 
tions of motion. For coalescing binary systems, the wave- 
forms of the emitted gravitational waves are expected to 
carry a great deal of information, and full exploitation 
of the expected observations will require accurate theo- 
retical models of the waveforms pT|. This requirement 
has prompted the computation of point-particle equa- 
tions of motion (as well as radiation reaction effects) to 
higher and higher post-Newtonian orders. Most recently 
the coalescence waveform's phase has been computed up 
to post-3. 5-Newtonian order see also Refs. |13t. il-^. 
At post-5-Newtonian order and higher, the concept of 
point-particle equations of motion will break down due to 
effects related to the finite sizes of the bodies J^] . How- 
ever an argument due to Damour jl^ indicates that the 
point-particle equations should be well defined at lower 
orders, up to and including post-4. 5-Newtonian. 

A third type of equation of motion applies to systems 
of N interacting bodies whose sizes cannot be neglected. 
These equations consist of the point particle equations 
of motion supplemented by tidal interaction terms. In 
principle, if one included tidal interactions to all multi- 
pole orders, and in addition coupled the equations of mo- 
tion to a dynamical description of the internal degrees of 
freedom in each body, one would obtain a complete de- 
scription of the system, equivalent to that provided by 
the continuum equations of motion (up to radiative ef- 
fects). 

For a system of bodies of typical size ^ i?, of typi- 
cal mass ~ M, and with typical separations ~ D, the 
force F that acts on one of the bodies can be written 
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schematically as^ 
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Here we use geometric units with G = c = 1. The terms 
inside the curly brackets are as follows. On the first line, 
the 1 is the usual Newtonian force between two point 
particles and the second and third terms are the post-1- 
Newtonian and post-2-Newtonian point-particle correc- 
tions. On the second line, the first term is the correction 
due to Newtonian tidal couplings; the minimum value of 
I allowed isl ~ 2 corresponding to quadrupolar coupling. 
The second term describes the post-l-Newtonian tidal 
couplings. Here the minimum allowed value I is lower 
than in the Newtonian case due to gravitomagnetic in- 
teractions which have no Newtonian analogs. This min- 
imum value is I = 1/2, corresponding to spin-orbit cou- 
plings (assuming that the bodies internal velocities are 
maximal, v ~ M/R). 

The purpose of this paper is to compute in detail the 
post-l-Newtonian tidal interaction terms in Eq. (|1.1|) . for 
all values of I, for a system of N bodies. The explicit 
form of these terms has not been derived before, although 
there is a substantial literature on this topic j2|, [l^ IitI 

El E 13, m m m m m 113 . 

There are a number of motivations for this computa- 
tion. First, as described in Refs. in the area of 
celestial mechanics future experiments and observations 
in the solar system will provide very high precision data. 
For example, there are current plans to increase the ac- 
curacy of lunar laser ranging from the current centimeter 
level to the millimeter level 28]. The future astrometric 
missions SIM (Space Interferometry Mission) and GAIA 
(Global Astrometric Interferometer for Astrophysics) are 
expected to measure angles to an accuracy of a few mi- 
croarcseconds, as compared to the current accuracy of 
milliarcseconds. In the radio, VLBI (Very Long Base- 
line Interferometry) observations currently can yield pre- 
cisions of order 10 microarcseconds (27|. Also, the pro- 
posed future laser astrometric test of relativity (LATOR) 
mission '29] would be sensitive to post-2-Newtonian ef- 
fects, and therefore would likely require detailed model- 
ing of post-l-Newtonian tidal effects. 

Second, gravitational wave measurements of coalescing 
binary compact stars will likely have some ability to de- 
tect finite size effects for sufficiently strong signals (SOj. 
Although post-Newtonian tidal effects will in many cases 
be small compared to Newtonian tidal effects, there are 
some situations where the post-l-Newtonian effects dom- 
inate. An example is the gravitomagnetic resonant exci- 



tation of Rossby modes in neutron stars that are spinning 
at ^ 100 Hz, which could be detectable with LIGO for 
moderately strong detected inspirals ■ 



B. Tidal coupling in post-Newtonian theory 

__The textbook treatment of post-l-Newtonian gravity 
4] is inadequate for the treatment of tidal interactions 
for several reasons, as explained by Damour et. al. 0]- 
First, the standard treatment uses a single global coor- 
dinate system. Although one can write down the contin- 
uum equations of motion for a given body in that coordi- 
nate system, it is very difficult to separate out the gravi- 
tational influences of the other bodies from the self-field 
of the body, since the fractional distortions of the coor- 
dinate system produced by the other bodies can be large 
even when the tidal distortion of the star is negligible. 
The development of approximation schemes such as lin- 
ear perturbations about an equilibrium state is hindered 
by the fact that the equilibrium state is not described in 
the usual way in the global coordinates. 

This difficulty has been comprehensively addressed 
in a series of papers by Brumberg and Kopeikin (BK) 
and by Damour, Soffel and Xu (DSX) 
0, lifl nil • These authors developed a detailed the- 
ory of post-l-Newtonian reference frames, in which each 
body has associated with it a coordinate system natu- 
rally adapted to that body. DSX also developed a for- 
malism to compute translational equations of motion in- 
cluding the coupling to all the mass and current mul- 
tipole moments of each body^. They applied their for- 
malism to compute equations of motion including spin 
and quadrupole couplings. In this paper, we extend the 
DSX results in two ways. First, by simplifying their for- 
malism we are able to compute the explicit form of the 
translational equations of motion, including all the multi- 
pole couplings. Second, we give a derivation that is valid 
for strongly self-gravitating objects as well as weakly self- 
gravitating objects^. We need only assume that the post- 
l-Newtonian field equations are satisfied in a weak field 
region surrounding each body. The bodies' internal gravi- 
tational fields can be arbitrarily strong; in particular our 
assumptions do not exclude black holes. By contrast, 
DSX assumed the global validity of the post-l-Newtonian 
continuum field equations, and so their derivation ap- 
plies only to weakly self-gravitating objects. Our result 
also generalizes existing derivations of the Newtonian 
and post-l-Newtonian |^ equations of motion for 
strongly self-gravitating objects that incorporate only a 



^ These scalings apply to generic bodies; if the bodies are spheri- 
cally symmetric the scalings are of course altered. 



^ The BK and DSX formalisms have recently been generalized to 

the parameterized post-Newtonian frame work for scalar-tensor 

theories of gravity by Kopeikin and Vlasov [33 . 
^ That is, we show that the dominant fraction al erro rs scale as 

0(M^/_D2); global post-Newtonian methods show only 

that these errors are 0{M'^ / B?) or smaller. 
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FIG. 1: An illustration of our assumptions for a system of 
bodies. Each body is surrounded by a strong field region 
which is excluded from our analysis. Surrounding these strong 
field regions are weak field buffer regions. Each body's center 
of mass worldline and mass and current multipole moments 
are defined in terms of the behavior of the metric in that 
body's bufi'er region. We assume that the vacuum post-1- 
Newtonian field equations are satisfied in all the buffer regions 
and in the regions of space between the buffer regions. 



few low order multipoles. Similar derivations to higher 
post-Newtonian orders including monopole terms only 
can be found in Refs. [sHllSB^. 

One of the key ideas underlying our derivation is that 
the equations of motion are determined entirely by the 
local field equations in weak field regions between the 
bodies. This was originally pointed out by Weyl and by 
Einstein and Grommer; see Thorne and Hartle and 
references therein. Each body is surrounded by a vac- 
uum, weak field region called a "buffer region" [SJI, and 
the quantities entering into the equations of motion are 
defined in terms of the behavior of the metric in those 
buffer regions (see Fig. QJ. In particular, our multipole 
moments are defined in terms of the behavior of the met- 
ric in the buffer regions. Our definition of multipole mo- 
ments is thus more general than the definition in terms of 
integrals over sources used by DSX. However, our multi- 
pole moments do coincide with those of DSX in the case 
of weakly self-gravitating bodies. 

Our derivation consists of a number of steps: (i) The 
definition of each body's current and mass multipole mo- 
ments and center-of-mass worldline in terms of the behav- 
ior of the metric in that bodies' buffer region, (ii) The 
definition for each body of a set of gravitoelectric and 
gravitomagnetic tidal moments that act on that body, 
again in terms of the behavior of the metric in that bod- 
ies' buffer region. For the special case of weakly self- 
gravitating bodies, our definitions of these multipole and 
tidal moments agree with definitions given previously by 
DSX. (iii) The derivation of a formula, for any given 
body, of the second time derivative of its mass dipole 
moment in terms of its other multipole and tidal mo- 
ments and their time derivatives. This formula was ob- 
tained previously by DSX for weakly self-gravitating bod- 



ies, (iv) A derivation of the relation between the tidal 
moments acting on each body and the multipole mo- 
ments and center-of-mass worldlines of all the other bod- 
ies. A formalism to compute this relation was developed 
by DSX; we simplify their formalism and compute the 
relation explicitly, (v) The deduction from the previous 
steps of the explicit translational equations of motion, 
whose form has not been previously derived. 

C. Results for equations of motion 

We next describe our results for the equations of mo- 
tion. We label each body by an integer A, with I < A < 
N. We use a harmonic coordinate system {t,x^) that 
covers all of spacetime except for the strong field regions 
near each body. The position of body A in this coordi- 
nate system is parameterized by a function = '""zf'{t) 
called the "center of mass worldline" . This function is 
defined precisely in Sec. IV CI below. It does not cor- 
respond to an actual worldline in spacetime; rather it 
parameterizes the location of the local asymptotic rest 
frame (see below) attached to the Ath body. That is, it 
is encoded in the behavior of the metric in a weak field 
region surrounding body A in the same way that the ac- 
tual center-of-mass worldline of a weakly self-gravitating 
body would be encoded. 

Associated with each body A is a coordinate system 
{satUa) which is defined only in that body's buffer re- 
gion, and which is adapted to the body in the sense that 
it minimizes the coordinate effects of the external gravi- 
tational field due to the other bodies as much as possible. 
This coordinate system is discussed in detail in Secs. lIII Dl 
and IV AI below. We will call the corresponding reference 
frame the "body frame" or, following Thorne and Hartle 
[S^ l- the body's "local asymptotic rest frame". The de- 
tails of the transformation between the body-adapted co- 
ordinates (sajUa) ^^'^ the global coordinates {t,x^) are 
important for the purpose of deriving the translational 
equations of motion. However, for the purpose of using 
the equations of motion, one only needs to know the fol- 
lowing. First, the time coordinate sa is a "proper time" 
associated with body A. It corresponds to the proper 
time that would be measured by an observer in the local 
asymptotic rest frame of body A. In that local asymp- 
totic rest frame it is related to the global frame time 
coordinate t by 

SA^SAit), (1.2) 

where the function SAit) is determined by a differential 
equation [Eqs. I|1.7d(l and H1.8(l below]. Second, the lead- 
ing order relation between the spatial coordinates y\ and 
x' is just a translation together with a time-dependent 
rotation [cf. Eqs. (|On|l and lf05|l below]: 

x'^-zfit) + uf^it)y^^. (1.3) 

The rotation matrix U^\t) describes dragging of iner- 
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tial frames^; a differential equation for its evolution is 
given below. The body- adapted coordinates {sA,y\) ro- 
tate with respect to distant stars, while the global coor- 
dinates (<, a;*) do not. 

Each body A has associated with it a unique set of 
mass multipole moments 



and 



(1.4) 



for I = 0,1,2... which are symmetric, tracefree, spa- 
tial tensors with I indices, of which mass dipole M^{sa) 
vanishes identically. It also has a unique set of current 
multipole moments 



StisA) = St,...a.{sA), 



(1.5) 



for I — 1,2,.... These quantities are functions of the 
body's proper time sa- In the absence of interactions 
with other bodies the mass monopole M"^ and the spin 
Sf are conserved. 

We will obtain below coupled equations of motion for 
the center of mass worldlines of all the bodies. Appearing 
in these equations as unknowns will be the mass multi- 
pole moments M^{sa) for I > 2, and the current mul- 
tipole moments S'^{sa) for / > 2. In order to obtain a 
closed system of equations, one would need to supplement 
the equations of this paper with equations determining 
the evolution of these multipole moments. We discuss 
further below various circumstances and approximations 
in which the evolution of the multipoles can be computed. 

Next, we define the moments M."^(t) and S^{t) to be 
the body's mass and current multipole moments, trans- 
formed to the non-rotating frame, and expressed as func- 
tions of the global time t. These moments are given by 
the equations [cf. Eqs. H5.35|) and (|5.36|l below] 

Mt..a,it) = UtHt)...Ut'mil„,,[sAit)], 

(1.6a) 

St..a,it) - L/.1'^'Ht)...C/a1'^'(054...,;[.A(t)].(l.6b) 

We can now write down the schematic form of the 
equations of motion. They can be written as 

<B k'aB k-.B c-B cSi 
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(1.7a) 
(1.7b) 
(1.7c) 

(1.7d) 



* The time derivative of this rotation matrix is actually of post-1- 
Newtonian order, so to Newtonian order this is a constant ma- 
trix. In the body of the paper we assumed that this constant, 
Newtonian order rotation matrix is the unit matrix. The descrip- 
tion of our results given here allows this constant matrix to be 
arbitrary; this slight generalization would be useful to describe 
systems that evolve for a time long enough that the accumulated 
rotation due to frame dragging is of order unity. 



= H 



(1.7e) 



Here Tf^^ T\ Tf, T and are functions of their 
argument whose specific forms are discussed below. In 
these equations the dependencies on the time derivatives 
Alf , 7W£ and Sf^ only occur for I > 2. Also the mass 
dipoles Mf vanish identically. Therefore, if we assume 
that the moments Alf{sA) and S^{sa) are known for 
I > 2, Eqs. (|l.t)a|) - (|1.7e|l form a closed set of evolution 
equations which can be solved to obtain the center of 
mass worldlines as well as the rotation matrices U"^ and 
time functions SA{t)- 

We remark that the three equations (|1.7a|l - H1.7c(l 
by themselves form a closed set of evolution equations 
for the variables '""zf{t),A4^{t), and Sf-{t), if we as- 
sume that the moments A4f^{t) and Sf{t) are known for 
I > 2. However, approximation schemes for computing 
the multipole moments for I > 2 usually yield the vari- 
ables Mf{sA), S^{sa) rather than the variables M^{t), 
Sf{t). This is because the moments Mf{sA) and Sf{sA) 
are the physical moments that would be measured by an 
observer in the the local asymptotic rest frame of body 
A. In such cases we must enlarge the set of variables 
°"^zf{t),A4'^(t),Sf{t) to include the rotation matrices 
U'^(t) and time functions s^(i) in order to obtain a closed 
set of equations. 

We also note that it is formally consistent to post-1- 
Newtonian accuracy to replace Eq. 
pier relation 5f (i) = S^{t). 



|1.6b|l with the sim- 
Nevertheless it can be use- 
ful in some circumstances to use the more accurate re- 
lation Ijl.Bbp . for example for systems which evolve for 
sufficiently long times that the rotation matrices U^"' 
become significantly different from unity. 

We now discuss the functions Tf, F^, Tf-, and 
Tf- that appear in Eqs. p.7a|l - (|1.7e|l . The functional 
form of is one of our key results. It is given by Eq. 
()f).ll(l below, with the coefhcients modified according to 
the substitutions given in Eq. (|6.19|l and in Appendix 
IFl These modified coefficients are obtained by combin- 
ing the results of this paper with those of the second 
paper in this series 37], which we will call paper II. The 
functions and Tf" are standard functions that can be 
derived from Newtonian stress-energy conservation for 
weakly self-gravitating bodies, and their explicit func- 
tional forms are respectively given in section IV of paper 
II 3 7] and in Eq. (jEH) below. The validity of Eqs. ifLTbl) 
and p.7c(l for strongly self-gravitating bodies is derived 
in paper II 37]. The function is given by [cf. Eqs. 
inSU, (E!2SaI), and below] 
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Here K is the multi-index bi . . . 6fc, vba = \ — 

„BA _ f cai^B cin^A\ /„ „_,j ^BA _ „BA ^BA 

n — [ z - z j/rsA, and — n^^ •■•"f,fc ■ 

Finally, the function is given by Eq. 16.20|l below. 

We next discuss various approximation schemes in 
which the equations of motion (|1.7a() ~ (|1.7e|l can be sup- 
plemented by methods for obtaining the evolution of the 
mass and current multipole moments M^{t) and Sf^{t) 
for Z > 2 in order to obtain a complete, closed set of 
equations. Some examples of such approximations are as 
follows. 

(i) The simplest case is when the effect of all the I > 2 
multipoles is negligible, and one can set Aif^ = Sf = 
for alH > 2. This yields the monopole-spin truncated 
equations of motion discussed in Refs. J,6, 34]. 

(ii) Another simple case is when the evolution of the 
multipoles of each body is dominated by dynamics inter- 
nal to that body, and is negligibly influenced by the tidal 
fields of the other bodies. In this case, one can solve for 
the evolution of the multipoles Adf^SA), Sf{sA) of each 
body separately, and then insert those multipoles into 
the equations of motion (|l.()a|) - 1)1. 7e|) . This application 
will be valid only if the timescale over which the bod- 
ies' multipoles evolve is sufficiently long [sll; see Sec. II Dl 
below for further discussion of this point. 

(iii) Another useful case to consider is that of rigid 
bodies. As noted by Thorne and Giirsel JS^], in general 
relativity a body's rotation can be rigid only if its angular 
velocity (with respect to its local asymptotic rest frame) 
is constant. If the angular velocity is changing, for ex- 
ample due to precession, then the body cannot be rigid 
due to Lorentz contraction effects. However, to linear or- 
der in the body's angular velocity the motion is rigid |39j| . 
The analysis of Thorne and Giirsel can be adapted to the 
present context, if the bodies' rotations are slow enough 
that they can be idealized as rigid. In this case, the time 
dependence of the mass multipole moments Mf{sA) for 
I > 2 can be parameterized in terms of a time-dependent 
rotation matrix U^'^{sa)- 

Mt..aMA) = Ut:^{sA)...Utr{^A)M-t..,.^. 

Here the moments are constant; these are the 

moments in the co-rotating frame which rotates with the 
body. We define the angular velocity r2^(syi) in the usual 
way as U^^U^'' = eacb^"^- Then, the co-rotating frame 
spin = U^°'S^ is related to the co-rotating frame 
angular velocity = U^^'^a via 39, 4^ 

Si{sA)^ltj,^i{sA). 

where is the (constant) moment of inertia tensor^. 
Similarly the higher order current multipole moments are 



^ Thorne and Giirsel |3fll have shown that for fully relativistic 
stars, as for Newtonian stars, the moment of inertia tensor is 
constant, independent of the angular velocity, up to linear order 
in the angular velocity. 



given by 

where I"^ _ is a higher order generalization of the mo- 

ai...aib ° » 

ment of inertia tensor [22|- Combining these relations 
with the equations of motion l|1.6a|l - (|1.7e|l yields a 
closed system of equations which can be solved for the 
center of mass worldlines, the rotation U^°'{sa) of each 
body with respect to its local asymptotic rest frame 
[sA, Ua), and the rotation U'^(i) of that local asymptotic 
rest frame with respect to distant stars. These equations 
describe torqued precession of relativistic objects, gener- 
alizing the free precession equations of Thorne and Giirsel 

(iv) For weakly self-gravitating bodies one can use the 
formalism developed by DSX ^ J6, 17] to obtain a post- 
1-Newtonian description of the internal dynamics of each 
body, for example by using post-l-Newtonian stellar per- 
turbation theory. Coupling such a description to the 
equations of motion yields a closed system of equations. 

(v) Lastly, for fully relativistic, spherical stars, one can 
compute the leading order effects of tidal interactions by 
combining the results of this paper with linear relativistic 
stellar perturbation theory using matched as ymp totic ex- 
pansions; see, for example, Refs. [iH li^. li^ l44j . For ex- 
ample, if one is interested only in the mass quadrupoles, 
and one restricts attention to the dominant, fundamental 
I = 2 modes with no radial nodes, then one has a relation 
of the form 

M^jisA) = J ds'^KisA - s'j,)Gt,{s'j,). 

Here K{sa — s^) is a Green's function which can be 
computed from stellar perturbation theory, and G^- is 
the body-frame gravitoelectric tidal moment that acts 
on body A, which is defined in Sec. IV Hl below and which 
can be computed in terms of the worldlines and multi- 
pole moments of the other bodies. Combining this rela- 
tion with the equations of motion (|1.7a() - H1.7c() again 
yields a closed system of equations, if one neglects the 
mass multipoles for I > 3 and the current multipoles. 



D. Domain of validity of our results 

As mentioned above, the key assumption which we 
make in deriving our results is that the post-l-Newtonian 
vacuum field equations are satisfied in a weak field region 
between the bodies; see Sec. IV Al below for more details. 
We are unable to give a derivation of this assumption 
from first principles. However, in this subsection we dis- 
cuss various physical effects which can cause our assump- 
tion to break down, and we make estimates of the sizes 
of these effects. We believe that the assumption should 
be generally valid aside from the effects discussed in this 
subsection. 
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The first type of correction are post-2-Newtonian cor- 
rections to the metric in the weak field, vacuum region 
between the bodies. These will give rise to fractional cor- 
rections of order M^/_D^, where M is a typical mass and 
D a typical separation of the bodies, cf. the third term 
on the first line of Eq. We can estimate as follows 

when these corrections will be larger than the tidal cou- 
pling terms which we retain. The estimate given in the 
last term of Eq. Hl.l|l can be refined by multiplying it by 
the dimensionless measure = A^i/(Mi?') of the Zth 
mass multipole. Demanding that this quantity be larger 
than the post-2-Newtonian, point particle term in Eq. 
yields the criterion 

^^^r{^)^R- (1-9) 

Thus, the post-2-Newtonian terms will always dominate 
at sufficiently large separations D, but for any given 
I > 2 there will be a range of values of D for which 
the post-l-Newtonian tidal terms dominate, as long as 
M/R. In particular this will be true for generic 
(non-symmetric) bodies for which e; ~ 1. Similar es- 
timates apply to current multipole couplings for I > 2. 
Thus, there is a nonempty regime in which the post-l- 
Newtonian tidal couplings computed here dominate over 
post-2-Newtonian, point-particle effects. 

Note, however, that this range of values of D gets 
smaller as the strength of internal gravity ~ M/i? in- 
creases. In the limit of Af ~ i? of a black hole, the post- 
2-Newtonian terms are always comparable to or larger 
than the post-l-Newtonian tidal terms. Therefore, our 
results cannot be applied consistently to black holes with- 
out including post-2-Newtonian and higher terms in the 
equations of motion. 

Another type of correction, which is also formally of 
post-2-Newtonian order, is that due to the time depen- 
dence of the mass and current multipole moments of the 
individual bodies 38]. The post-l-Newtonian solutions 
[Eqs. (|3.5a|l ~ (|3.5c|l below] do not exhibit the correct re- 
tarded dependence on these moments, that is, they are 
functions of Mi(t) and SL{t) rather than Afi(i — r) and 
SL(t — r). If these moments vary on a timescale r, then 
the corresponding fractional corrections to the mass mo- 
ments scale as D'^/t'^, and the fractional corrections 
to the current moments scale as D'^/t^. Demand- 
ing that these corrections be smaller than the post-l- 
Newtonian accuracies of these quantities M/D and 
~ 1 respectively) yields the criterion r ^ 13 for the cur- 
rent moments, and the more stringent criterion 

r»D[-) (1.10) 

for the mass moments. Fractional corrections to the post- 
l-Newtonian tidal interactions will be of order unity for 
r ^ D{D /MY/^ . The criterion t ^ D essentially says 
that all of the bodies lie in the near zone of the gravita- 
tional radiation produced by any one body, and not in 



the wave zone p8] , and the criterion (|1.10() is a somewhat 
stronger requirement than this. 

To illustrate the criterion (|1.10() it is useful to consider 
some examples. First, if the time evolution of a body's 
moments is driven by tidal interactions with other bod- 
ies, then the timescale is of order r ~ D(Z3/M)^/^, and 
the criterion is satisfied. Second, suppose that we have 
a 3-body system consisting of two black holes in a tight 
binary, together with a third body in orbit around the 
binary. We can model such a system as a 2 body system 
using the formalism of this paper, treating the black hole 
binary as a single body^ whose mass and current mul- 
tipole moments are evolving with time due to internal 
dynamics. Then, early in the gravitational- wave driven 
inspiral of the black hole binary, the criterion (|1.1U|) will 
be satisfied and our results for the equation of motion 
will be valid. As the black hole binary gets tighter how- 
ever, eventually the orbital period will become shorter 
than H1.10|l and our results will no longer be applicable. 

This second example illustrates that the post-l- 
Newtonian approximation can sometimes completely 
break down, even in the supposed weak-field region be- 
tween the bodies. During the final coalescence of the 
black holes the gravitational radiation metric perturba- 
tion will become temporarily as large as the Newtonian 
potential in the region between the binary and the loosely 
bound companion. Our results are not applicable to such 
systems, in which one of the bodies emits a strong burst 
of gravitational radiation. Further work is required to 
deduce the form of the translational equations of motion 
in this type of situation. 

There are two other assumptions made in our deriva- 
tion which slightly restrict the domain of validity of our 
results. First, we assume that a coordinate system which 
covers the weak field region between the bodies can be 
smoothly extended to cover the bodies' interiors (see Sec. 
nil Bl below) . This assumption essentially restricts the 
spatial topology of the bodies' interiors, and excludes 
objects like eternal black holes, wormhole mouths and 
naked singularities. It does not exclude realistic, astro- 

tsical black holes for the reason explained in Sec. lIII I^ 
. Second, in order for a body's multipole moments to 
be definable, it is necessary that there exist two concen- 
tric coordinate spheres surrounding the object, such that 
the region between the two spheres is vacuum, in a par- 
ticular coordinate system centered on the body (see Sec. 
nil Bl below). This assumption might break down when 
two bodies get within one or two radii of one another, 
slightly before they actually touch. 



^ Although our formalism cannot be consistently applied to indi- 
vidual black holes unless supplemented with post-2-Newtonian 
and higher order point particle terms, our formalism can be ap- 
plied to black hole binaries treated as a single body. This is 
because binaries are less compact than black holes. 
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E. Organization of this paper 

The structure of the paper is as foUows. In Sec. Ill Al we 
introduce our notations for the post-l-Newtonian contin- 
uum field equations, and following DSX we define a class 
of gauges (conformally Cartesian harmonic gauges) that 
we use throughout the paper. Section 111 ti\ presents a 
simplified version of the theory of post-l-Newtonian ref- 
erence systems of Rcfs. H El 113 IH IH H El • The 
key result of this section is the explicit parameterization 
(|2.17|l of the residual gauge freedom within conformally 
Cartesian harmonic gauge in terms of a number of freely 
specifiable functions of time and one harmonic function 
of time and space. 

Section IIIII is devoted to the definitions of the mass 
multipole moments Afi(t), current multipole moments 
SL{t), gravitoelectric tidal moments GL{t) and gravito- 
magnetic tidal moments HL{t) associated with a given 
object and a given conformally Cartesian, harmonic co- 
ordinate system. These definitions are given in Sec. lIII Bl 
in terms of the general solution l|3.5a|l - H3.5c|l of the 
post-l-Newtonian field equations in a vacuum region be- 
tween two concentric coordinate spheres that surround 
the object (the object's "buffer region"). Section Fill CI 
analyzes how all of these moments transform under the 
class of allowed gauge transformations discussed in Sec. 
HlBl In Sec. HITdI we describe gauge specializations that 
fix the gauge freedom completely and accordingly deter- 
mine the multipole and tidal moments uniquely. We call 
the resulting coordinate system a body-adapted coordi- 
nate system. Section IlII El gives a definition of multipole 
moments and tidal moments associated with a given ob- 
ject, a given worldline and a given coordinate system. 
These moments arise only in intermediate steps in the 
derivations of this paper and not in our final results. Fi- 
nally, in Sec. IIII Fl we compare the moment definitions 
used here with other definitions in the literature. 

Section Hvl derives the law of motion for a single body, 
that is, the relation between the second time derivative of 
its mass dipole moment and its other multipole and tidal 
moments and their time derivatives. The assumptions 
and result are described in Sec. IIV Al A general descrip- 
tion of the surface-integral method of derivation which 
we use is given in Sec. IIV Bl In Sec. lIVCl we give some of 
the post-2-Newtonian vacuum field equations which are 
needed for the derivation. Section llV DI derives the single 
body law of motion to Newtonian order, as a warm-up ex- 
ercise. Finally, the post-l-Newtonian derivation is given 
in Sec. IIV El This derivation uses an idea due to Thorne 
and Hartle ,34j to deduce the value of a complicated sur- 
face integral from previous weak-field computations of 

DSX 0,131. 

Section lays the foundations for treating a system 
of N interacting, finite-sized bodies. Our assumptions 
are described and discussed in Sec. IV Al In Sec. IVBI we 
define, for each body, a set of body-frame multipole and 
tidal moments associated with that bodies' adapted co- 
ordinate system. These are the moments that would be 



measured by an observer in that bodies local asymptotic 
rest frame. Section IV CI defines the configuration vari- 
ables that specify the location, orientation etc. within 
the global coordinate system of the local asymptotic rest 
frame which is attached to that body. These variables in- 
clude the center of mass worldline and also the time func- 
tions and rotation matrices discussed in Sec. II CI above. 
In Sec. IVDi we define for each body multipole and tidal 
moments associated with the global coordinate system. 
These quantities appear only in intermediate steps in our 
computations and not in our final results. The relation 
between the global-frame multipole and tidal moments 
and the body-frame multipole and tidal moments is com- 
puted in Sec. IV El Section IV Fl defines the modified ver- 
sions Ml and Sl of the body-frame multipole moments, 
discussed in Sec. II CI above, which are defined with re- 
spect to a frame that is non-rotating with respect to dis- 
tant stars, and which are expressed as functions of the 
global time coordinate. These are the moments that ap- 
pear in the final equations of motion. 

Finally, Sec. I VII gives the derivation of the complete, 
explicit translational equations of motion for the N body 
system. 



F. Notations and conventions 

Throughout this paper we use geometric units in which 
G = c = 1. We use the sign conventions of Misner, 
Thorne and Wheeler ; in particular we use the metric 
signature (— , -|-, -I-). Greek indices (/i,z/ etc..) run 
from to 3 and denote spacetime indices, while Roman 
indices (a, b, etc..) run from 1 to 3 and denote spatial 
indices. The spacetime coordinates will generically be 
denoted by (a;°,a;*) = [t.x"^). Spatial indices are raised 
and lowered using 5ij, and repeated spatial indices are 
contracted regardless of whether they are covariant or 
contravariant indices. We denote by the unit vector 

/r, where r = = ^ bijx'^xK 

When dealing with sequences of spatial indices, we use 
the multi-index notation introduced by Thorne [3^ as 
modified slightly by Damour, Soffel and Xu 01 . We use 
L to denote the sequence of I indices a\a2 ■ ■ ■ ai, so that 
for any Z-index tensor^ T we have 

Tl = Taiaa-.-aj- (l-H) 

If Z = 0, it is understood that Tl is a scalar. If Z < 
then Tl = 0. Wc define L — 1 to be the sequence of / — 1 
indices aia2 . . . a;_i, so that 



Here by "tensor" we mean an object which transforms as a ten- 
sor under the symmetry group SO(3) of the zeroth order spatial 
metric Sij , not a spacetime tensor. 
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If Z = 0, then by convention T^-i = 0. We also define N 
to be the sequence of n indices aia2 ■■■<!„, and L — N to 
be the sequence of / — n indices a„+ia„+2 ■ ■ - ai, so that 
we can write a relation like 



(1.13) 



as Gl — SnTl-n, for any tensors G, S and T. We de- 
fine K, P and Q to be the sequences of spatial indices 
6162 • ■ • ^fc, C1C2 . . . Cp, and did2 ■ ■ ■ dg, respectively. Re- 
peated multi-indices are subject to the Einstein summa- 
tion convention, as in SlTl- We also use the notations 



and 



dL = da^c 



(1.14) 



(1.15) 



We use angular brackets to denote the operation of 
taking the symmetric trace-free (STF) part of a tensor. 
Thus for any tensor T^, we define 



T<i> = STFL(ri). 



(1.16) 



where STF^ means taking the symmetric trace- free pro- 
jection on the indices L. For example, li I — 2, we have 

See Appendix ^ for the general definition of STF pro- 
jection, and for a collection of useful relations involving 
STF tensors. 

Throughout this paper, symbols will generally denote 
functions (as is common in mathematics) rather than 
physical quantities (as is common in physics). For ex- 
ample, in Sec. we define a mass multipole moment 
M^{sa) which is a function of a time coordinate sa- 
In that section we also use a different time coordinate 
t. Then, the notation Mf{t) will always mean AIl{sa) 
evaluated at Sa — t, rather than M^[sa{t)]. 

Finally, for the aid of the reader an index of symbols 
is provided in Table UTI 



II. POST-l-NEWTONIAN CONTINUUM FIELD 
EQUATIONS AND GAUGE FREEDOM 

In this section we summarize the form of the post-1- 
Newtonian field equations that we use, and analyze the 
residual gauge freedom left after imposing our gauge con- 
ditions. Our notation follows closely that of Weinberg 
[47| ■ thou gh we relate our conventions to those of DSX 

Ulillillii. 



A. Metric expansion and field equations 

In the post-l-Newtonian approximation, one considers 
a one-parameter family of solutions to Einstein's equa- 
tions of the form 

ds^ = -[l + 2e^<^> + 2e'^{<i>^ +ip)+0{e^)]{dt/e)^ 



+ [2e^Q+0{e^)]dx'{dt/e) 



-Ik 



0{e'^)]dx'dx'. 



(2.1) 



Here the Newtonian potential <I>, the post-Newtonian po- 
tential -0, the gravitomagnetic potential Q and the spatial 
metric perturbation hij are functions of the coordinates 



^0 



t and a;*, but are independent of the parameter e. 
The corresponding expansion of the stress-energy tensor 
is 



(2.2) 



where ■T^'' is the Newtonian-order piece and p'T^'' is 
the post-l-Newtonian order piece. The post-Newtonian 
expansion parameter e used here is equivalent to the ex- 
pansion parameter used in some other treatments 
(we use units in which G = c = 1). 

We note that many presentations of the post-l- 
Newtonian equations use a time coordinate t that is re- 
lated to our time coordinate t by 



t 



(2.3) 



In (t, x') coordinate systems the various powers of e that 
appear in the expansions (|2.1(l of the metric and H2.2|l 
of the stress energy tensor differ from those given here. 
The (t, x') coordinate systems have the advantage that 
the pointwise limit as e — *■ of the metric exists and is 
a flat, Minkowski metric, but have the disadvantage that 
the potentials $, V' and Q become e-dependent. That e- 
dependence is usually accounted for by inserting an extra 
factor of s whenever a time derivative of a potential is 
taken. By contrast, in the coordinates used here, the 
potentials are independent of e and no such extra factors 
of e are needed. 

Assuming the validity of Einstein's equations for the 
metric (|2.1|l and stress energy tensor (|2.2|) for all e in 
some open interval < e < Eq then implies that the 
metric 6ij + e^{hij + 2^6ij) is flat to O(e^) Therefore 
one can always choose coordinate systems^ in which hij = 
—2^Sij, for which metric expansion simplifies to 



ds^ = 



-[l + 2£2$ + 2e4($2+^) 
+ [2e^Q + 0{e'^)]dx\dt/e) 
+ [Sij - 2£2$^y + 0{e'^)]dx'dx\ 



0{e')]{dt/ef 



(2.4) 



* This statement is always true locally, and is true globally if the 
spatial domain of the coordinates is simply connected. In this 
paper the spatial coordinate domains will always be simply con- 
nected. 
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Such gauges are called conformally Cartesian 0; we 
will restrict attention to conformally Cartesian gauges 
throughout this paper. 

We will also assume the harmonic gauge condition 

df^iV^gn = 0, (2.5) 

which at post-l-Newtonian order and for conformally 
Cartesian gauges reduces to 



^9$ 



dC 



= 0. 



(2.6) 



Harmonic coordinate systems are usually conformally 
Cartesian, but one can have local coordinate patches 
which are harmonic but not conformally Cartesian. 

The metric (|2.4f) . stress-energy tensor H2.2f) and gauge 
condition (|2.6|l 

imply the standard^ harmonic-gauge^*^ field equations 

El 



4^ TOO 
167r 



47r '■■T" 



(2.8a) 
(2.8b) 
(2.8c) 



For most of this paper we will be concerned with the 
vacuum versions of the field equations H2.8a|l - (|2.8c|l . 

These are 



0, 
0. 



(2.9a) 
(2.9b) 
(2.9c) 



For later use we note the expansions of the electric and 
magnetic components £ij = i?oiOj and Bij = —^eikiRkWj 
of the curvature tensor. In vacuum regions these tensors 
are symmetric and traceless and can be expanded as 



-2 pnC 
O7. 



and 



(2.10) 
(2.11) 



^ If we add £2 times Eq. |Z8^ to Eq. IT^ we can combine these 
two field equations into the single wave equation —d^w/dP + 



-47rcr -I- 0(e ), where w ■ 



-<I> — e^i/) and 



(2.7) 



This is the notation used by DSX The potential w satisfies a 
fiat-space wave equation, with hidden e dependencies that come 
from the definition of w. For our purposes it will be more useful 
to keep the expansion in e fully explicit by using the two elliptic 
equations l^.ijal and I2.8bl for <I> and ip instead of the above 
single hyperbolic equation for w. 

In the "standard post-Newtonian gauge" |^ the gauge condition 
12.61 is replaced by S^" -I- V ■ = 0, and the field equations 
iZ8bl and IZSct are replaced by V^V = ^■K^"r^° + iirT" and 
V2(» = levrTO' -I- <i-,i. 



Here "fy = is the Newtonian electric tidal tensor, 
and 



<,»^j> + 4$$,<jj> -I- (2.12) 



is the post-Newtonian electric tidal tensor. The angular 
brackets denote an STF projection, cf. Sec. II Fl above. 
The post-Newtonian magnetic tidal tensor is 



(*>i)' 



where 



B = V X C- 



(2.13) 



(2.14) 



is the so-called gravitomagnetic field. 

For later use, we also note the definition of the gravi- 
toelectric field E used by DSX |2|: 



(2.15) 



B. 



Parameterization of residual gauge freedom in 
conformally Cartesian harmonic gauge 



In many applications of post-l-Newtonian theory, the 
elliptic harmonic-gauge field equations l|2.8a|) - H2.8c|) are 
valid in all of space. In such cases one normally solves 
the field equations by imposing the boundary condition 
that all the potentials go to zero at spatial infinity. This 
boundary condition determines a unique solution to the 
field equations and a unique choice of gauge. 

However, in this paper we will be dealing with situa- 
tions where the field equations H2.8a|l ~ (|2.8c() do not have 
unique solutions, due to residual gauge freedom. There 
are two reasons for this residual gauge freedom ,2]. First, 
even in cases where the field equations are valid in all of 
space, there is in fact no physical reason for imposing 
that the potentials go to zero at spatial infinity. Instead, 
the physical boundary condition to impose is that the 
components (I2.1U|) and 12.11|1 of the Ricmann curvature 
tensor go to zero at spatial infinity. One then finds that 
there is a large class of solutions of the harmonic gauge 
field equations, when the sources are fixed. This non- 
uniqueness is present even at Newtonian order; there 
are many solutions to the Newtonian Poisson equation 
(|2.8al) with the boundary condition ^ as r ^ 00. 
In the Newtonian context, the new solutions are simply 
the standard solution transformed to accelerated frames. 
Similarly, in the post-Newtonian context, the additional 
solutions correspond to the original solution transformed 
to reference frames that are accelerated, rotating or oth- 
erwise modified with respect to the standard reference 
frame. 

In Sec. nil Dl below, when considering a system of N 
interacting bodies, we will need to construct a coordinate 
system adapted to each body. Exploiting the additional 
freedom of allowing accelerated, rotating coordinate sys- 
tems will be crucial for our construction of those adapted 
coordinate systems. 
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The second reason for residual gauge freedom in the sit- 
uations considered in this paper is that we wiU be consid- 
ering spacetimes containing strong-field regions in which 
the post-Newtonian approximation is not valid. There- 
fore, we must analyze the field equations (|2.8a() - H2.8c() 
on some spatial region T) which is not all of space. In this 
case, the boundary conditions imposed on the potentials 
<I>, ip and Ci on the boundary dT> of T) influence the solu- 
tion. Part of the information inherent in those boundary 
conditions is gauge, and part of the information is phys- 
ical. 

DSX 0, Kopeikin ^ and Klioner and Voinov 
have derived a complete parameterization of the residual 
gauge freedom present in some region T) of space after 
the conformally Cartesian condition has been imposed; 
see Ref. |20] for a review. In this subsection we give a 
simplified, streamlined version of the DSX analysis in our 
somewhat different notation. We also specialize the DSX 
analysis by imposing in addition the harmonic gauge con- 
dition. 

Our starting assumptions are as follows. We assume 
the existence of two different coordinate systems (i, x*) 
and (f, x*) on V, each of which is conformally Cartesian 
and harmonic. We also assume that both coordinate sys- 
tems are such that the metric admits an expansion of the 
form (|2.4II . In particular, there exist potentials x^), 
-ipit.x^) and Q{i,x^) such that 



-[l + 2e2|. + 2£'^($2+^) 
-[2e^Q + 0{e'^)]dx\dt/e) 
-[(5y - 2e^m.,j + 0{e^)]dx'dx^ 



Oie')Kdt/s) 



(2.16) 



In Appendix ^ we show that the most general relation 
between the two coordinate systems that is compatible 
with these assumptions is"'^^ 

= x' + z\t)+e^h\i,x^) + Oie'^), 
t = t + e^a{t,x^) + e^f3{t,x^) + 0{e^), (2.17) 



and where overdots mean derivatives with respect to the 
time argument. 

We next discuss the meaning of the various freely speci- 
fiable functions adt), z^{t), h\{t), Rk{t), and f]h{t,x^) 
that appear in the coordinate transformation (|2.17|l - 
(|2.18cll . At Newtonian order there appears the function 
of time ac{t), which governs the normalization of the time 
coordinate at 0{£'^). In standard treatments of Newto- 
nian gravity, this freedom is fixed by the usual assump- 
tion $ ^ as r ^ oo. In the present context, however, 
this coordinate freedom is not fixed. There also appears 
the spatial 3- vector z*(i), which parameterizes the trans- 
lational motion of the new frame with respect to the orig- 
inal frame, to Newtonian order. At post-l-Newtonian or- 
der, one has the 3-vector hl{t), which parameterizes the 
post-Newtonian translational motion of the new frame 
with respect to the original frame. There also appears the 
spatial 3-vector Rk{t), whose time derivative is an angu- 
lar velocity that parameterizes the slow, post-Newtonian 
rotation of the coordinate system. Finally, there is the 
function Phit, x^) which governs the normalization of the 
time coordinate at 0{e*). This function is not completely 
freely specifiable but must be harmonic, i.e. 



V2/?h = 0. 



(2.19) 



It is straightforward to compute how the potentials 
transform by combining the metric expansions l|2.4|l and 
(|2.1(j|) with the coordinate transformation given by Eqs. 
(|2.17l) - (|2.18c|l and using the tensor transformation law 
for the components of the metric. The results are 



- Zi z- 



y3 



$(i,X^) 

Cz(f, x') = x') - [4$(i, x') + a]^^ 



dhj 
' di 



dx^ 



dp 
dx^ 



(2.20a) 



(2.20b) 



and 



where 

a{t, X-') 
h'{t,x^) 



ac{t) + Xiz\t), 



(2.18a) 



(2.18b) 
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^{t)xk + -ac{t) 



Mt,x^), 
(2.18c) 



Up to constant displacements in time and time-independent spa- 
tial rotations. We also assume that the coordinate transforma- 
tion is orientation-preserving and time-orientation preserving. 



tpit^x^) = ■il!{t,x^) - Q{t,x^)z'' + a 



d<P{t,x^ 



+2<!>it,x')z,z' - {ai' - h') 



dt 

d^(t,x^) 
dx^ 



1 / . , . \ 2 . . , . . , dhi 8(3 1 . 2 
+"^^^-^^ + W-2"- 



Here the function $(t, x^) is defined as 
$(i,x^') = $[t,x^' -l-z^'(t)]. 



(2.20c) 



(2.21) 



where the right hand side is evaluated at the point 
= -f z^{t), t = i. We define the functions ipit,^:^) 
and C'(t, a;-') similarly. 

The transformation laws H2.20a| - H2.20c|l are ex- 
pressed in terms of the functions a(t, x^), P(t,x^) and 
h\t,x^) defined by Eqs. ifTTTjl - (|2.18c|l . More explicit 
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versions of the transformation laws, in which they are ex- 
pressed in terms of the freely specifiable functions ac{t), 
z^{t), hl.it), Rkit) and /3h(f, S-'), can be obtained by sub- 
stituting the definitions (|^T7|) - 12.18cl) into Eqs. H2.20a|) 
- ^HM . 

There are two special subgroups of the group (|2.17|) 
of transformations that will be of importance later. The 
first subgroup applies when the spatial domain T) is all of 
space, and when in addition one imposes the boundary 
condition that all the potentials vanish at spatial infinity. 
In this case it is easy to show that only uniform relative 
motion of the two frames is allowed, = ft,^ = 0, that 
/3h = 0, Qc = and that is constant. These well- 

known "post-Galilean" transformations are discussed in, 
for example. Sec. 39.9 of Ref. 

The second important subgroup is the subgroup pa- 
rameterized by the harmonic function /3h, for which 
Qfc = = /i* = = 0. The corresponding coordinate 
transformations are 

x'=x' + 0(e^) , t = f+e4/3h + 0(£^), (2.22) 

and the potentials transform according to 

One can show that the Newtonian and post-Newtonian 
pieces of the connection coefficients F^^ and of the Rie- 

mann tensor components R^iSj^ are invariant under this 
subgroup. As noted by DSX 0, this subgroup corre- 
sponds to a gauge freedom in post-l-Newtonian theory 
analogous to that of electromagnetism^^. 



III. DEFINITIONS OF AN OBJECT'S MASS 
AND CURRENT MULTIPOLE MOMENTS AND 
TIDAL MOMENTS 

A. Overview 

As discussed in the introduction, a crucial part of the 
derivation of the equations of motion for strongly self- 
gravitating bodies is the definition of an object's mass 
and current multipole moments and also tidal moments 
in terms of the behavior of the metric in a weak field 
region surrounding the object. In this section we discuss 
the definitions of these quantities. 

For orientation, we start by reviewing the definition 
and status of multipole and tidal moments in Newtonian 



If one requires only that the coordinate systems be conformally 
Cartesian and not harmonic, then the most general coordinate 
transformation is still given by Eqs. STT7\ - I2.18cl . but with the 
modification that the function /3h can be arbitrary rather than 
being harmonic 0. 



gravity. Suppose that in some reference frame {tjX"^), 
there exists a region T) between two spheres of the form 

r_ < |a;| < r+, 

for some radii r_ and r+, in which there are no sources. 
Then the Newtonian potential $ satisfies the Laplace 
equation (|2.9a|) in V. The general solution for $ in D 
can then be written in terms of a multipole expansion as 

°° f-iV+i 1 1 

/=o ■ II- 

(3.1) 

Here L denotes the multi-index ai 02 ... a; and de- 
notes the tensor x^^x^^ . . .x'^' , cf. Sec. II I' I above. The 
quantity "M^ (i) is the Zth Newtonian mass multipole mo- 
ment associated with the reference frame {t, x^ ) of the ob- 
ject or objects in the region \x\ < r_. [The superscript 
"n" in "Mi denotes "Newtonian" .] Similarly the quan- 
tity X^Lit) is the Ith Newtonian tidal moment associated 
with the reference frame {t,x^) that acts on the region 
r < due to sources outside r — r^, where r = \x\. 
The moments "Ml and "Gl are both STF tensors. 

The expansion (|3.1|l can be taken as the definition of 
the moments "Afi and °Gl; it is possible to invert Eq. 
(|3.1|) to obtain explicit expressions for these moments in 
terms of surface integrals of $ in the domain V (see Ap- 
pendix^. If we additionally assume that the Newtonian 
Poisson equation l|2.8a|l is valid everywhere in r < r_, 
then we obtain the conventional formula for the mass 
multipole moments as an integral of the Newtonian mass 
density T"°: 

"ML{t)= f T™(i,x^)x<^>d3x. (3.2) 

Here the angular brackets denote an STF projection, cf. 
Sec. II Fl above. As is well-known, the field-based defini- 
tion (|3.1() of the multipole moments is of greater gener- 
ality than the integral-based definition 13.2|l . since the 
former is applicable to strong field sources that possess 
an asymptotic region in which the Newtonian description 
is a good approximation. 

Next, we recall that the moments °Mi(i) and "GL(t) 
depend on the choice of reference frame or coordinate 
system (t, x^). They change when one switches from one 
reference frame to another according to x x + z{t) 
[cf. Eq. (|2.17l) above, to Newtonian order]. This ambi- 
guity is conventionally resolved in Newtonian physics by 
specializing to the reference frame in which the origin of 
coordinates coincides with the center of mass of the ob- 
ject, or, equivalently, in which the mass dipole moment 
"Mi(i) vanishes. 

Consider now the corresponding situation at post- 
l-Newtonian order. A suitable generalization of the 
integral-based definition H3.2|l of multipole moments for 
an isolated system was given by Blanchet and Damour 
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[48j . That definition was generalized to the case of sev- 
eral interacting bodies by DSX DSX also gave a 
field-based definition of multipole moments analogous to 
the definition Ij^.H) of "Ml, and gave a different type 
of definition of tidal moments. In this section we will 
review, simplify^'^ and generalize the definitions of DSX 
. Our analysis will be more general than theirs because 
we will consider spacetimes where the post-l-Newtonian 
field equations are not satisfied everywhere, whereas DSX 
assumed the global validity of those field equations. 

As is the case at Newtonian order, the post-Newtonian 
multipole and tidal moments associated with a body are 
not uniquely defined but depend on the choice of refer- 
ence frame or coordinate system. This ambiguity can be 
resolved, as in Newtonian theory, by making a specific 
choice of canonical reference frame adapted to a given 
body. However, the freedom in choice of reference frame 
is much larger at post-Newtonian order than at Newto- 
nian order, cf. the discussion in Sec. Ill Bl above. There- 
fore the specialization to a body-adapted frame is more 
involved. 

The remainder of this section is organized as follows. 
In Sec, nil Bl we give a form of multipolar expansion of the 
general solution of the post-l-Newtonian field equations 
that serves to define the multipole and tidal moments 
associated with a given body and with a given coordi- 
nate system. Section Fill CI discusses the gauge transfor- 
mation properties of the moments, and Sec. IIIIDl defines 
the body-adapted gauge that fixes the moments uniquely. 
Finally in Sec. lIIIEl we generalize our definition of multi- 
pole moments to define moments associated with a given 
coordinate system about a given specified worldline. 

B. Definition of mass and current multipole 
moments and tidal moments 

To define the multipole moments of a body, we start 
by assuming the existence of a local coordinate sy;stem 

13 The DSX definitions, given in Eqs. (6.9a) - (6.10b) of Ref. Q. in- 
volves a splitting of all the post-Newtonian potentials into pieces 
associated with the individual objects. Our simplified version of 



{t,x^) with following properties, (i) The range of the 
coordinates contains the product of the open ball 

\x\ < r+, (3.3) 

where r-|_ is some radius, with some open interval {to, ti) 
of time, (ii) The vacuum post-l-Newtonian field equa- 
tions H2.9a|l - (|2.9c(l are valid in a spatial region V of the 
form 

r_ < \x\ < r+, (3.4) 

for some non-zero radius r_. 

These assumptions allow us to define the multipole and 
tidal moments of a body or bodies in the region r < r_ , 
where r = |a;|. In parallel with the Newtonian case dis- 
cussed above, the second assumption allows for the possi- 
bility of strong field sources for which the post-Newtonian 
approximation is not valid in the region r < r_. When 
applying this definition to systems of several bodies, we 
will choose both r_ and r+ to be of order the distance 
between the bodies. 

The first assumption, that the coordinates in the region 
V can be extended into the interior to cover the body, 
is not actually necessary for the definition of multipole 
moments. However it will be used later in the derivation 
of equations of motion so we include it here. Note that 
this first assumption does not exclude the possibility that 
one or more black holes could reside in the region r < 
r_. As pointed out by Thorne one merely needs to 
choose as the t ~ constant surfaces an appropriate set of 
time slices that pass though the interior of the collapsing 
object(s) that form the black hole(s). 

Given these assumptions, the general solution of the 
vacuum field equations (|2.9a|) - (|2.9c|l in the region V 
can be expanded in terms of STF tensors as 

their definitions do not require any such splitting. 



°° (-1)'+1 1 1 



1=0 



E 

1=0 

1 



(-1) 



i+i 



"Mi(t) 



(2/ + 1) 



1 \x\^ 



Al(0 



1 , (-1)'+^ \x\ 



-[^T;.(t)-..(i)].---^^'t;.(0.n, 



E 

1=0 

1 

~n 



(-1) 



'+1 r 4 . 4/ 21-1 



4(2?- 1) 



L{t) + Y^—^eji^aiHL-l>]{t) ^y— — "G'<L_i(t)(5a,> 



(3.5a) 



(3.5b) 



(3.5c) 
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The quantities that appear in these equations are the 
following. First, there are the Newtonian mass multipole 
moments "ML{t) and tidal moments "Ghit) that were 
discussed in Sec. lIII Al Second, there are post-Newtonian 
corrections ^"ML{t) and ""GLit) to these quantities. [The 
superscripts "pn" denote post-Newtonian.] We shall call 
the quantities 



ML{t) = "ML{t) 



and 



GUt) = XJL{t)+e^'"GL{t) 



(3.6) 



(3.7) 



the total mass multipole and tidal moments, respectively. 
Third, there are current multipole moments SL{t), and 
a new set of tidal moments (t) related to gravitomag- 
netic forces. Following DSX, we will refer to Hj^it) and 
Ghit) as the gravitomagnetic and gravitoelectric tidal 
moments, respectively. Fourth, there are moments fJ-L{t) 
and I'Lit) that contain information about the coordi- 
nate system being used, but do not contain any gauge- 
invariant information about the body. We shall show 
below that it is always possible to find a gauge in which 
liL{t) = VL{t) — 0. We shall call these quantities "gauge 
moments" . 

All of the quantities "A/l, "Gl, ^"A/l, ''"Gl, Ml, vl, 
Sl and Hl are STF on all their indices. As in the New- 
tonian case, the expansions (|3.5a() - (|3.5c() can be taken 
as the definition of all of these moments for a given coor- 
dinate system; one can invert these expansions to obtain 
explicit expressions for all the moments in terms of sur- 
face integrals of various combinations of derivatives of 
the potentials [cf. Appendix IE] . The moments "Ml, 
"Gl, ""A/l, ^'Gl, fJ-L are defined for all I > 0, whil e z/l, 
Sl and Hl are defined only for Z > 1. In Eq. (|3.5b|l and 
throughout this paper, it is understood that 1/^=0 for 
I = 0. 

We shall call the pieces of the potentials that would 
diverge as \x\ — > oo the tidal pieces of the potentials. 
Correspondingly, we will use the phrase "tidal moments" 
to refer to any of the moments "Gl, ""Gl, Hl and 
that appear in the coefRcients of the growing terms in 
Eqs. (|3.5a() - (|3.5c() . These moments encode the gravita- 
tional influence of other bodies on the body in the region 
I a; I < We shall call the remaining pieces of the po- 
tentials, which involve the moments ''Ml, ^"Ml, Sl and 
^L, the intrinsic pieces. 

We now turn to the derivation of the expansions (|3.5bp 
and H3.5c|l . We start by writing the general solution in 
the region T) of the Laplace equation H2.9c|) for the grav- 
itomagnetic potential as 

°° (-l)'+^„ .... 1 1- 



1=0 



(3.8) 



Here the quantities ZiL and I^l are STF on their L in- 
dices only, and not on the i index. Next, we insert the 



expressions (|3.8ll and H3.5a|l for the gravitomagnetic and 
Newtonian potentials into the harmonic gauge condition 
()2.6|1 . This gives the relations^"* 



"Af = 0, 



J<iL> — 



l + l 



and 



-An J 



(3.9) 
(3.10) 

(3.11) 



Next, we use the identity given in Eq. (6.21) of Ref. 
that allows one to express in terms of irreducible, STF 
tensors any tensor that is STF on all its indices except 
its first index. Specifically, if T^l is any tensor satisfying 
TiL = Ti<^L>, then we have 



T.L - TI + '^ + ^,^<a;^^\y, + '5.<a,Tni, (3-12) 



.(0) 



.(-1) 



where 



(3.13a) 



(0) _ 



1 + 



J^Tjk<L~l^ai>jk (3.13b) 



and 



T 



(-1) _ 



L-l 



21-1^ 
21 + 1' 



(3.13c) 



In order to apply this identity to the tensors ZiL and I^l , 
we define 



Sl — —^Zjk<:L-i^ai>jk, 

Hl = Yjk<L-l<^ai>jk, 

VL = y<L> 



and 



Ml = Z. 



(3.14a) 

(3.14b) 
(3.14c) 

(3.14d) 



We now insert the relations (|3.10() and (|3.11() and the def- 
initions H3.14a|l - (|3.14d(l into the general identity (|3.12|l 
to obtain 

7 - 4 • 41 ^2?-l 

(3.15) 



Thus, the time independence of the Newtonian mass monopole 
"M can be derived either from stress-energy conservation at 
Newtonian order in the interior of the body (for a weakly self- 
gravitating system), or from the validity of the post-Newtonian 
vacuum field equations and harmonic gauge condition in the far 
field of the body. This type of phenomenon, where one can avoid 
dealing with the interior physics by going to one higher post- 
Newtonian order in the far field, is well known in the literature 
on equations of motion, and will be encountered again in Sec. 

ITVl 
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and 



I 



l + l^i^<-' — ^^^ 21 + 1 

(3.16) 

Inserting these formulae into the expansion 1)3. 8|l yields 
Eq. ^X^ . 

We remark that the parameterization (|3.8(l of the grav- 
itomagnetic potential in terms of the tensors I^l and ZiL 
will frequently be more convenient to use in our computa- 
tions below than the fully STF parameterization (|3.5c|l . 

Consider now the post-Newtonian potential ip. We can 
write the general solution in the region "D of the vacuum 
field equation Ij2.9bp as the sum i/j = i/jp -I- i/jh of a particu- 
lar solution ipp of the inhomogeneous equation and a gen- 
eral solution i/ih of the homogeneous equation V^V — 0- 
A particular solution can be obtained by inspection, us- 
ing the expansion l|3.5a|l of the Newtonian potential: 



HL-l>j - 



4(2/ - 1) ^ 



°G<L-lSa. 



l>i- 



1=0 ^ 



2 /!2(2/-h3) 



(3.17) 



The homogeneous solution can be written in a form par- 
alleling the expansion H3.5a|l of the Newtonian potential: 



/=0 

Next we define 



^^^-"A/lSl^-tt^-Glx^. (3.18) 



and 



(/ + l)(2/ + 3)' 



(3.19) 



(3.20) 



Inserting the definitions (|3.19l) and (|3.20|l into the ho- 
mogeneous solution (|3.18|) and adding the particular so- 
lution ifTTTjl yields Eq. (|33H)l . The reason for choosing 
the particular parameterization given by Eqs. H3.19|l and 
is so that the moments ''"^/^(t) and ""GLit) be 
invariant^^ under the subclass 12.22|l of gauge transfor- 
mations; see Sec, nil (Jl below for more details. 

We now specialize to the situation, considered by DSX, 
where the post-l-Newtonian field equations with sources 
(|2.8a|) - H2.8c|l are assumed to hold for all r < r+, i.e. 
in the interior of the body. In this special case, we now 
show that our definitions of the quantities ^"M^, ^"Gl, 
Hl and 5*^ are equivalent to the DSX definitions of these 
quantities. 

Consider first the multipole moments p°Afi and 5*^. 
Equations (6.9a) and (6.9b) of DSX "^l define the mo- 
ments Mi(i) and ^^(t) in terms of an expansion of "lo- 
cally generated" pieces of the potentials. Their splitting 



of the potentials into "locally generated" and "external" 
pieces is defined by their Eq. (4.5), and is easily seen to 
be equivalent to the splitting which we discussed above 
of our expansions l|3.5a|l - (|3.5c|l into intrinsic terms and 
tidal terms. Therefore it is sufficient to show that the 
DSX expansions (6.9a) and (6.9b) coincide with the tidal 
terms in our Eqs. (|3.5a|) - H3.5c|) . This follows from the 
definition H3.6|) and the relations W — — $ — e^ip and 
W, = -Ci/4 between the DSX potentials {W,W') and 
our potentials tp and C*- 

Consider next the gravitoelectric and gravitomagnetic 
tidal moments Gl and Hl. These are defined by DSX in 
terms of STF projections of gradients of the external (or 
tidal) pieces of the gravitoelectric and gravitomagnetic 
fields evaluated at the origin of spatial coordinates, cf. 
Eq. (6.13) of Ref. 2]. Inserting our expansions H3.5a|l - 
(|3.5cll into the definitions (|2.14() and (|2.15() of the grav- 
itomagnetic and gravitoelectric fields yields for the tidal 
pieces (denoted by a superscript T) of these quantities 



oo 



1=0 



HiLX 



<L> 



and 

eT 



(3.21) 



oo . 

Eyy i'G^L+e^'"G^L)^ 

1=0 



.<L> 



71-4 



2(2/4-3) 



21 + 1 



I + 1 



(3.22) 



These expressions agree to 0(£^) with corresponding ex- 
pansions derived by DSX [Eqs. (6.23) of Ref. 0], when 
the definition (|3.7|l is used. It follows that the two defi- 
nitions of tidal moments are equivalent. 

Thus, our definition of both the multipole moments 
"Ml, ""Ml, Sl and the tidal moments "Gl, ""Gl and 
Hl in terms of the general solution H3.5a|l - (|3.5c|l of 
the post-l-Newtonian vacuum field equations unifies and 
simplifies the two different types of definition given by 
DSX. Our definitions also generalize the DSX definitions 
to strong field sources. 

Finally, we note that in the case considered by DSX, 
one can alternatively define the post-Newtonian multi- 
pole moments by integrals over the source that are anal- 
ogous to the Newtonian integral (|3.2() Translating 
Eqs. (6.11) of Ref. 2] into our notation gives for these 
integrals 

Sl= f e^''<'''x^~^>^ T^'^ d^x (3.23) 



and 
""Ml 



3^3 



X-' X- 



00 



Except for ""Gl for I = which is not invariant. 



2(2/ + 3) 
(/-H)(2/ + 3) ^ '""^ i"' 



24) 
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Gauge transformation properties of the 
moments 



harmonic function /?h as 



In this section we compute how the various moments 
transform under the general transformation (|2.17|) from 
an original coordinate system (t,x^) and a new coordi- 
nate system {i,x^). Under that transformation, the spa- 
tial domain V defined by Eq. H3.4I) is mapped onto the 
domain 



r- < \x 



m < 



(3.25) 



to zeroth order in e. We therefore restrict the set of 
coordinate transformations to those that satisfy 



max \z{t)\ < 
t 



(3.26) 



This restriction ensures that the image V of the domain 
T) contains a non-empty region of the form f_ < \x\ < f+ 
for some radii r_, f+ with r_ > 0, and thus allows us 
to define multipole and tidal moments in the new coordi- 
nate system. We also parameterize the freely specifiable 



/3h 



1=0 



x\ II 



TL{i)2 



arbitrary smooth function 



in V 



in 



(3.27) 

where Pint is the image of the domain Q < r < r^. Here 
the quantities Xhit) and Ti(f) are STF on all their in- 
dices. This choice guarantees that the coordinate trans- 
formation (|2.17|l is defined and smooth on the entire do- 
main < r < r+, while maintaining the harmonic prop- 
erty of the coordinates in the domain < r < r+ where 
the post-l-Newtonian equations are valid. 

In the barred coordinate system we define the trans- 
formed moments °ML(f), "Ghit), ""MLit), ""GLit), 
%L{t), Z,L{t), SL{t), HL{t), VL{i) and JiL{i) by the fol- 
lowing barred versions 
and lEHSIl: 



J 



and 



E 

i=0 



1 1 



'ML(i)ai— --"GL(i) 



E 

i=0 



(-1) 



l+l 



\x\ l\ 

(21 + 1) 
(/ + l)(2/-h3) 



/tL(f) 



112(21 + 3) 



E 

4 



1 1. 



7.<L> 



l + l 



"M.Lit) 



Al - 21-1 _ 



l+l 



21 + 1 



4(2?- 1) 



YzLit) = l^iLit) + y-^ejj<a,-ffL-l>i(t) - ^^^^ 

I 



<L-1 



(3.28a) 



(3.28b) 
(3.28c) 
(3.28d) 

(3.28e) 



As before dots denote derivatives with respect to the 
time argument. By substituting the multipole expan- 
sions H3.5a|l — (|3.5c|l of the original potentials into the 
transformation formulae H2.20a|l — (|2.2()c|l . and compar- 
ing the results with the multipole expansions (|3.28a|) — 
, we can derive transformation laws for the various 
moments. 



To illustrate this procedure, we begin with the trans- 



formation of the Newtonian potential. We obtain 



^i,x'') = ^t,x'') + z,{t)x'--mz\t)+d,it) 



(=0 



x + z{t)\ 



'\xiL{t)[x + z(i)]^ -I- z,(t)x' - \z\i)k{^- 

(3.29) 
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Next we use the expansion 

1 °° 1 1 

4- r ril 



' ' n— 



X 



(3.30) 



which is vahd for \x\ > \z\ [cf. Eq. (|3.26|) above], and we 
compare with the muhipole expansion 13.28a|) . This gives 
the well known transformation laws for the monopole and 
dipole 



"M = "M, 



and 



"Mi = "Mi - "Mz„ 
together with the general result 



"Mr. - E(-i: 



71 + 1 



II 



n=0 



n\{l-ny. 



"M^nzl 



<nzl-n>- 



(3.31) 
(3.32) 

(3.33) 

+2 • • ■_£ai> J 



Here "M^nzl-n> denotes "M^ai...a„Za„^j^Za^ _ 
cf. Sec. II F\ above. Similarly we obtain from Eq. 
the transformation law for the Newtonian gravitoelectric 
tidal moments 



oo 



fe=0 



Here 



Af = 



(3.34) 



(3.35a) 
(3.35b) 



and Af = for ? > 1. Following DSX 0, we caU the 
quantities Af inertial moments. 



We similarly compute the transformation laws for the 
multipole and tidal moments Zn. and . We insert the 
expansion l|3.8|l into the transformation law l)2.20b|l and 
use the formulae H2.18al) - (|2.18c|l for the quantities a, 
and (3. We also use the parameterization H3.27|l of the 
harmonic function /3h. We then equate the result to the 
multipole expansion H3.28c(l and obtain the transforma- 
tion laws 



I 



ZiL — 



(3.36a) 



n\[l — ny. 

n=0 ^ ' 

-4Zi"M^NZL~N>) + l5i<ai^L-l>, 
oo ^ 

(3.36b) 



k=0 



Here 



Af 



AS 



ZiZj Z^iZj'^ "4^ Cijf^R}^ ^~ '2Sij Zj^Z}^ 



(3.37a) 
(3.37b) 

6 

Ay-fe = --Si<j'zk>, (3.37c) 
and all the other gravitomagnetic inertial moments A^^ 
are zero. 

The transformation laws for the fully STF moments 
Sl, Hl, and vl parameterizing the gravitomagnetic 
potential [cf. Eq. (|3.5c|l above] are obtained by splitting 
the expressions ()3.36all and Ij3.36b|l for ZiL and YiL into 
their fully STF pieces using the identities (|3.12|) . HA5p . 
l18II and lAlQI). The results are 



Sl = E(-l)" 



+1 



n=0 

n 



n\{l-ny. 



n{l + 1) [l-n) ■ 

!^<NZL-N> — ,f , ^•^ Zs"Mr<NZ(L-l)-N^ai>rs 



l{n + l) 



l{n + l) 



+ jZr ''Ms<N-lZ(L~l)-{N-l)^ai>rs H J ZrZg "-M<ArZ(i_ i) _7veQ, >rs 



(3.38a) 



Mi 



Ef + 1) 

^ ^ nUl + l-nV. 

n=0 ^ ' 



l+l-n ( 2 
' 2n - 1 - 



; + l (? + l)(2n+l) 



2/ + 1 



[i-<N-\ZL-{N-\)> 



4(^ + 1 -n)(2Z-2n + l) ■ 4n(/ + 1 - n)(2/ + 3) 

" (/ + l)(2/ + l)(n + l) '^r^h<NZL-f,> - + X){n + + 1) ''^^v<N-^^i^-^)-(.^-^)^-^>v3 



4n 



4G + 1-^). . 



■y^jTY^j "A^j<A'-iZL_(Ar-i)> ^■'^^ "M<nzl-n> 

8n{l + l-n) . 4:{l + 1 - n){l - n) 

-Zj Mj^N-lZ(L-l)-{N-l)Za,> H _^ ;^^^2; + 1) ^^^^ M<NZ(L-l)-NZai> 



il + l)i2l + l) 

4(/ + 1 - n){l - n) 



(? + l)(2^ + 3) 
21 + 1 



(3.38b) 
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and 



^ k\ 

k=0 



H 



lkz 



K 



8k ■ 



(3.38c) 



k=0 



+ 



4fc(2/ + 2fc- 1) 
{I + k){2l + 2k + I) 
Akl 

{I + k){2l + 2k + 1) 



Z<i'GL>K-lZ^ ^ + 



l + k + 1 

Ak{k- 1) 



[l + k){2l + 2k + I) 



Z'' £jm<iHL>mK-lZ^ 



'^K-l<iL-lZai>Z 



K-1 



r 



(3.38d) 



It is also possible to derive from the expressions H3.36a|) 
and l|3.36b|) transformation laws for the moments "Mn, 
and "GiL that enter into the fully STF parameterizations 
(|3.15|) and H3.16|l of Z^^ and Yn^ . The resulting transfor- 
mation laws are consistent with the transformation laws 
(jT^ and derived earlier for "Ml and "Gl; this 



consistency is an important check of the formalism. 

Finally we turn to the moments ^"AI^ and ^"Gl pa- 
rameterizing the post-Newtonian potential ip. The trans- 
formation laws for these moments are by far the most 
tedious to compute. The results are 



where 



and 



-"Mr = 



(; + i)(2; + 3)' 



ri=0 



n\{l-ny. 



'^<nZl-n>, 



""Mr 



21 + 1 



(? + l)(2/ + 3) 



+ 



/2U 




\2l^ 





fiL - z.jZjL + '2zjZj "Ml + (ofc - zjZj) "Ml - lU^ai "Ml- 

21 



i> 



(/ + 2)z, "M,L - Uj, "Ml ~ I "M,<L-iC/a,>j + l^f-p^U.^a, 'Ml-i>j + z, "M, 



1 



--^5ijZk + Z[i5j]k, 



(3.39) 



(3.40a) 

(3.40b) 
(3.40c) 

(3.40d) 
(3.40e) 



where 



p'l 



i^L + Y, ^{PLK + lz<ajL^i>K)z'^ + Tl - UAf ^ 
k=0 



PL = "'Vl-VL- Z,Y,L+2ijZj"GL + Uj'G!jL+lU.j<a,'XiL~l 



>J 



- l)z<,, "Gl-i> + lz<a, XiL-i> + (ac - Zjz') "Gl + ^(^f^Pi 
"Gl + 2zj "GjL + ZjZk °Gl + Zj "GjL, 

li^J^j) ^2' 



. ■ 1- . • _ 3 . .. . _ 



(3.41) 



(3.42a) 
(3.42b) 
(3.42c) 

(3.42d) 
(3.42e) 
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with all the other being zero. 

Note that the expression (|3.39|l for the transformed 
moment ^"1'/^ depends on the transformed moment /i^. 
However, using the transformation law (|3.38bl) for /i^ we 
can write the expression H3.39|l entirely in terms of the un- 
transformed moments, and the dependence on /ii cancels 
out. Similarly the expression H3.41|l for the transformed 
moment '"Gi depends on the transformed moment P^, 
but that dependence can be eliminated using the trans- 
formation law l|3.38d|l for vi^. 

Finally we note that the left hand sides of Eqs. H3.33|l . 
(ESI, (ESnai, IISEQ, ll^O^ - and iTTTni are 

functions of the new time coordinate t. The right hand 
sides are expressed as functions of t by evaluating the 
untransformed (unbarred) moments, which are functions 
of at t = i. This replacement of t by t in the argu- 
ments of the untransformed moments is implicit in the 
transformation laws (|2.20a|l - H2.20c() for the potentials'^. 



D. Specialization to body-adapted gauge 



The construction in Sec. IIII 51 above defined the mass 
multipole moments "A/i(i), ^"ML{t), current multipole 
moments SL{t), gravitoelectric tidal moments 'iJL{t), 
^XjLit) and gravitomagnetic tidal moments HL{t) that 
are associated with a given body and with a given choice 
of coordinate system. In this section we discuss how to 
obtain a unique set of multipole and tidal moments asso- 
ciated with a given body by specializing to a coordinate 
system that is adapted to the body in a certain way. 
The particular coordinate system defined here is rela- 
tively well known; see for example Ref. 0- We shall call 
it the "body-adapted gauge" . 

We start by reviewing the well-known construction at 
Newtonian order. From the transformation laws (|3.32|) 
and (|3.34l) for "Af i and °Gl, we see that the freedom car- 
ried by the Newtonian-order worldline z{t) could be used 
in either of two ways. One can set to zero either the mass 
dipole moment "Miit) or the I — 1 gravitoelectric tidal 
moment Gi{t). The second choice is not very useful, as 
it makes the worldline of the origin of spatial coordinates 
follow the tidal field instead of following the body. The 
first choice is the conventional and useful choice; choosing 
Zi{t) = °M,{t)/°M achieves 



(3.43) 



and the coordinates are then mass-centered to Newtonian 
accuracy. Normally such mass-centering would mean 
that the origin of coordinates coincides with the body's 



Note in particular that we do not evaluate the untransformed 
moments at t = t{t), where t{t) is the function t{t,x^) of the 
coordinate transformation I2.17i evaluated at = 0. The cor- 
responding correction terms have already been included in the 
derivation of Eqs. l2.2Uai - l2.2Ucl 



center of mass. In the present context, however, this con- 
clusion is not valid, since the origin of coordinates is in 
the strong-field region r < r_ of space where the New- 
tonian equations are not necessarily valid. Moreover the 
definition of the center-of-mass in the present context is 
somewhat subtle; we defer to Sec. IV CI the discussion of 
this definition. 

The second gauge specialization we make at Newtonian 
order is to set 



"G = 



using the choice 



adt) 



dt 



ZjZj + 



1=0 



(3.44) 



(3.45) 



cf. Eqs. (|3.34() and (|3.35a(l above. In Newtonian physics, 
this choice simply corresponds to adjusting the zero of 
the gravitational potential, which does not affect the dy- 
namics. 

We now turn to a discussion of the gauge specialization 
at post-Newtonian order. Without loss of generality we 
can assume that we have already achieved the Newtonian 
conditions "Mi = "G = 0, and so we can specialize to 
the purely post-Newtonian subgroup of the coordinate 
transformations which is characterized by (t) = adt) = 
0. That subgroup is parameterized by the functions hl(t), 
Rf^ (t) and by the STF tensors Al {t) and tl [t) that define 
the harmonic function /3h via Eq. H3.27|l . 

We first discuss the gravitoelectric sector. For the 
purely post-Newtonian coordinate transformations the 
transformation law H3.39|l for the post-Newtonian mass 
multipole simplifies considerably to 

-"Ml - ""Ml -1"M^l-iK^> ~ l'M,<L-iea,>jkR'' ■ 

(3.46) 

If we specialize to the mass dipole, the last term in Eq. 
(|3.46|) vanishes since the coordinates are already mass- 
centered to Newtonian order. This gives 



p°M, = -"Mi - "Mh^ 



(3.47) 



Therefore, as in the Newtonian case, we can use the 
translational freedom encoded in the function /ij to mass- 
center the coordinates by making 



-"Mi = 0. 



(3.48) 



The transformation law H3.41|l for the gravitoelectric tidal 
moments similarly simplifies to 



-"Gl= -"GL+hyV,L-V'Gj<L-iea,>jkR''+hi, (3.49) 

where the quantities are given by A'-^ = f, Kf = —hi, 
and A^ = for I > 2. By choosing r(t) suitably we can 
make 



-"G = 0, 



(3.50) 
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which is analogous to the Newtonian condition l|3.44() . 

Consider next the gravitomagnetic sector. For purely 
post-Newtonian coordinate transformations, the trans- 
formation laws H3.38a|) and H3.38c|l for the current mul- 
tipole moments Sl and gravitomagnetic tidal moments 
Hl simplify to 



Sl 
Hl 



Sl, 
Hl 



(3.51) 
(3.52) 



where Af = -2Ri and A|f = for Z > 2. Therefore we 
can choose to make 



(3.53) 



by an appropriate choice of the angular velocity R'^, as 
noted by DSX 0. This gauge specialization makes the 
I = part of the tidal piece of the gravitomagnetic 
field vanish. The resulting coordinate system slowly ro- 
tates relative to distant stars, in such a way that the 
leading order Coriolis acceleration due to the tidal grav- 
itomagnetic field is effaced. 

At this stage, the remaining coordinate freedom is pa- 
rameterized by the STF tensors Xl for Z > and tl 
for / > 1, which appear in the formula H3.27|l for the 
harmonic function /3h- Note that these tensors do not 
enter into the transformation laws H3.46|l . (|3.49|l . (|3.51(l 
and for the moments ^"Ml, ""Gl, Sl and Hl ^'^ ■ 

Therefore the values of these moments will be the same 
in all coordinate systems that satisfy the conditions spec- 
ified so far, Eqs. (ESJ, EiSl), and 
In other words, these moments (as well as the Newtonian 
moments "M^ and "Gl) are already uniquely defined by 
the conditions we have specified so far. Nevertheless, it 
is useful for some purposes to fix the remaining gauge 
freedom. To do this we consider the gauge moments \xl 
and vl- For purely post-Newtonian transformations the 
transformation laws (|3.38b(l and Ij3.38d|l for these mo- 
ments reduce to 



Mi 

VL 

where A J' — —h\ 
we can make 



(/ + l)(2? + 3) , 
P-L H „, , -, Al, 



Vl - TL 



2/ + 1 



I > 1, 



(3.54) 
(3.55) 



and A^ = for ? > 2. 



Vl 







It follows that 



(3.56) 



by choosing tl suitably, for alH > 1. Similarly we can 
make 



TABLE I: Free functions in the coordinate transformation 
and their role in defining the body-adapted harmonic gauge. 



Free function 



Role 



Kit) 

r{t) 
Rk{t) 
Mi), 1>1 
Al(F), />0 



sets "M, = 
sets "G = 
sets P"A/i = 
sets P"G = 
sets Hi — 
sets i>L = 
sets flL = 



ML =0 



(3.57) 



by choosing Xl suitably, for all Z > 0. We shall call 
the unique coordinate system-'^* that achieves all of the 
condi tions jXl^ . (|S3H1), EHOl, EiSSl), and 

(|3.57() the body-adapted coordinate system^^. 

To summarize, we have demonstrated that there is 
enough coordinate freedom to accomplish the following: 

1. Mass-center the coordinate system to post-1- 
Newtonian accuracy by setting "Mi = ^"Ali = 0. 

2. Set to zero the ^ = pieces "G, "t? of the tidal 
pieces of the potentials <I> and tp 

3. Set to zero the gravitomagnetic tidal moment Hi. 

4. Set to zero the all the gauge moments vl and fiL- 

The role of each free function in the coordinate transfor- 
mation H2.17|l in the derivation of the above conditions is 
recapitulated in table Q] 

We note that our definition of the body-frame coordi- 
nates or local asymptotic rest frame differs slightly from 
that of Thorne and Hartlc 34] . We require that the mass 
dipole Mi should vanish for all time, whereas Thorne and 
Hartle instead require the gravitoelectric tidal moment 
Gi should vanish for all time. At a given initial instant, 
they also demand that Mi and Mi vanish, while Mi can 
be nonvanishing. 

We next consider the special case treated by DSX, 
where the post-l-Newtonian field equations H2.8af) - 
(|2.8c|l are assumed to hold all the way down to r — 0. 
In this case the gauge freedom is somewhat reduced, 
since the function /3h must now be both harmonic and 
smooth for all r < r+. This requirement eliminates the 
terms parameterized by Al in Eq. (|3.27|) . as those terms 
diverge at r = 0. Therefore we no longer have sufficient 
gauge freedom to set to zero the gauge-moments /zl via 
Eq. lESlI). We still obtain a unique coordinate system 



Achieving this gauge invariance was the reason for picking the 
particular choices of parameterization of Eqs. 13.191 and I3.20i 
above. 



Up to constant displacements in time and time-independent spa- 
tial rotations, cf. Sec. Ill b1 above. 

This gauge is called the skeletonized-body harmonic gauge by 
DSX 13 . 
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by imposing the remaining requirements (|3.43() . H3.44(l . 
(E3H1), (|^3n|) . and but the ^z, moments 

wiU now in general be nonvanishing. This modified 
version of the body-adapted gauge will be important in 
Sec. HVI below. 



E. Definition of multipole and tidal moments 
about a given worldline 

In Sec. nil Bl above we defined the multipole and tidal 
moments of a body associated with a given coordinate 

By a worldline we mean simply a function = z^{t) which 
transforms appropriately under the group 12.171 of coordinate 
transformations. If z^{t) lies outside of the domain of definition of 
the coordinates then it does not correspond to an actual worldline 

I 



tpH, X 



where 



and 



The expansions (pnHa|l a nd |3.58b^ have the same form 
as Eqs. H3.5a|l and (|3.5b|l except that x is replaced ev- 
erywhere by X — z{t), and the various time derivative 
operators are allowed to act on the factors of x — z{t) as 
well as on the moments. Note that the formulae (|3.58d(l 
and H3.58e(l for ZiL and YiL contain extra terms involving 
compared to the original formulae (|3.15|) and H3.16|l . 

We now discuss the derivation of the expansions 
(|3.58a|l - H3.58cl) . First, one can verify that these expan- 
sions satisfy the field equations H2.9a|l - (|2.9c|l . Next, in- 
serting the expansions H3.58a|) and (|3.58c|l into the gauge 
condition H2.6|l yields the Newtonian conservation of mass 
equation (|3.9|l as before, and also the following replace- 



system. Those moments can be interpreted as being mo- 
ments about the origin a; = of that coordinate system. 
In this section, we generalize that definition to define 
tidal and multipole moments about a specified world- 
line^*^ a;* = and associated with a given coordinate 
system This more general definition will be used 

in Sec. EEl below. 

We assume that the vacuum post-l-Newtonian field 
equations are satisfied in a region of the form r_ < \x — 
z(t)\ < r+. Our definition of the moments is given by 
the following multipole expansions of the potentials $, 

and ijj in this region: 

in spacetime. See Sec. IV CI below for further discussion of this 
point. 



(3.58a) 



(3.58b) 
(3.58c) 



(3.58d) 



(3.58e) 

I 

ments for Eqs. H3.10|l and H3.11|l : 

Z<^L> = Y^"M^L + '^Z<^"ML> (3.59) 

and 

Y,jL = -4"GL-K4i,=GjL. (3.60) 

If we now define the moments Sl, Hl, and vl in 
terms of Zn, and Yn, using the same formulae (|3.14a() 
- (|3.14d|l as before and use the decomposition identity 
(|TT^ . we obtain Eqs. H3.58dp and H3.58e|) . Finally, the 
expression (j3.58b|l for -0 is chosen so that the moments 



1^0 

oo 



\x~z{t)\ II 

1 (2/ + 1) d 

\x-z{t)\ ^ {I + 1){21 + 3)di 



\x-z{t)\ 



(-iy+1 9^ 



\x~zit)\ 



1 



4|K(i)[--^(t)]^}-2,(2^^3)^,. 



{'•GUt)[x-zit)]^\x~z{t)\^}y 



°° ( lV+i 1 



1 



1=0 



\x-z{t)\ II 



Y,L[x-z{t)Y 



4 • Al 21—1 



4(2/ - 1) 



4(2/ + 1) 



2/ + 3 
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"Ml, p"Ml, "Gl, ""^l, Hl and S'i defin ed by these ex- 
pansions are invariant under the group H2.22|l of gauge 
transformations. This invariance can be verified by us- 
ing a parameterization of /3h of the form (|3.27|l with x 



replaced hy x — z{t). 

For later computations it is useful to expand the time 
derivatives in Eq. (|3.58b(l and express results in terms of 
STF tensors. This computation gives 



J 



where 



E 

/=0 



1 



X ~ z{t)\ 



1 



2/ + 3 



\x-z{t)\ 



21 + 1 



Fl[x - z{t)]^ + Jl 



\x-z{mx^z{t)Y 



2(2/ + 3) 



(3.61a) 



21- 1 
2/ + 1 



Zj "MjL + ZjZj "Ml + 2zj "Mjl + 2/i-,i<a, "Ml~i>j 



(3.61b) 



Pl = "Ml + 2/i<„, ■Ml-i> + l{l - l)i<a,ia,_i "Ml-2> + lz<a, °A^L-i>, 



(3.61c) 



and 



1 



21 + 3 



lz-'z-' 



21 + 1 



<L-lZai>, 



(3.61d) 



Jl 



L - 2 "G 



'VLjkz'z'' - X^LkzK 



r 



(3.61e) 



All of the tensors Nl, Pl, Fl and Jl are STF, while the 
tensors Zil and YiL are STF only on their last I indices, 
as before. 

Finally, we note that there is no natural, unique defini- 
tion of multipole and tidal moments about a given world- 
line associated with a given coordinate system. We have 
chosen a particular definition, but there are other defini- 
tions that are equally valid. For example, if we replace 
Eq. H3.58d|l with 



A ■ Al 

"M,L - -j-^£ji<aiSL-l>j + AZi "Mdio 



l + l 



21-1 ^ , , 

+ ^^_^ ^ 0»<a,ML-l> +4'Mj<L_iZa,>, (3.62) 



then field equations and harmonic gauge condition are 
still satisfied. However, while the moments defined by 
this equation are still invariant under the group (|2.22l) 
of gauge transformations, they differ from the moments 
defined by Eq. Ij3.58d|l . By contrast, the definitions of 
multipole and tidal moments about the origin of coordi- 
nates discussed in Sec. lIII 51 above are essentially unique. 
However, this lack of uniqueness will be unimportant for 
our purposes, since multipole and tidal moments about 



a given worldline will appear only in intermediate steps 
in our computations and not in our final results. 



Comparison with other definitions of multipole 
and tidal moments in the literature 



As we have discussed, our multipole moments coin- 
cide with those of Blanchet and Damour 48] and of DSX 
0, 0, 01 for weakly self-gravitating bodies. For iso- 
lated systems these moments agree to leading order with 
the unique moments that can be defined for stationary 
systems in general relativity :49], and with the asymp- 
totic radiative multipole moments of Thorne [3^ 1 , in the 
appropriate limits, as noted by DSX 

For non-isolated systems, there is another notion of 
mass and current multipole moments, defined in terms 
of the metric in a buffer region surrounding a body, due 
to Thorne and Hartle 0|. These moments are defined 
in full general relativity by the same type of surface in- 
tegrals as used here [cf. Appendix |E|, but applied to the 
full metric rather than to the post-Newtonian potentials. 
These moments depend on the choice of coordinate sys- 
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tem used to evaluate the surface integrals, but the mag- 
nitude of the resulting ambiguities can be estimated and 
in many applications are small enough to be unimportant 
[s^l- By contrast, the multipole moments used here are 
defined only the context of post-l-Newtonian theory, but 
are unique. 

Our tidal moments also coincide with those of DSX 
for weakly self-gravitating bodies. They also appear to 
coincide with the tidal moments defined by Suen in the 
context of stationary systems in full general relativity 




IV. POST-l-NEWTONIAN LAWS OF MOTION: 
A SINGLE BODY 



A. Overview 



of time, (ii) On the spatial region T> given by r_ < 
r < r-i-, for some r_ > 0, the coordinates are conformally 
Cartesian and harmonic. Also the vacuum Einstein equa- 
tions are valid on D for a one-parameter family of metrics 
of the form of Eq. I|4.11|l below. Essentially this says that 
the Newtonian, post-l-Newtonian and post-2-Newtonian 
vacuum field equations are valid on T). The reason for 
imposing the post-2-Newtonian field equations in addi- 
tion to the post-l-Newtonian field equations is discussed 
below. 

These assumptions allow us to define the multipole 
moments °Ml(<), p"ML(i), SL{t) which characterize the 
sources in the region r < r_ , as well as the tidal moments 
'^^(t), "'X^Lit) and HL{t), as discussed in the previous 
section. They also imply some formulae relating time 
derivatives of the multipole moments. At Newtonian or- 
der these formulae are 



As discussed in the introduction, the derivation of 
equations of motion for several interacting bodies can be 
divided into two pieces: (i) A derivation of a formula'^^, 
for any given body, of the second time derivative of its 
mass dipole moment in terms of its other multipole and 
tidal moments and their time derivatives, (ii) A deriva- 
tion of the relation between the tidal moments acting on 
each body and the multipole moments and center-of-mass 
worldlines of all the other bodies. In this section we will 
carry out the first of these tasks. The second task will 
be the subject of Sees. IVl and IVII below. 

We start by describing the assumptions and the result. 
As in Sec. IIII Bl above, we assume existence of a local 
coordinate system (t,x^) with following properties: (i) 
The range of the coordinates contains the product of the 
open ball 



\x\ < r+, 



(4.1) 



°M(<) ^ 0, (4.2a) 

OO ^ 

"M,(t) = ^-"Afi(t)"G,i(i). (4.2b) 



where is some radius, with some open interval (to, ti) 

DSX call this formula the "law of motion" 



1=0 



Here the first formula is just the Newtonian conservation 
of mass derived earlier [Eq. (|3.9|) above]. The second 
formula equates the acceleration of the center of mass to 
the acceleration produced by the external tidal moments 
coupling with the multipole moments of the system. This 
formula is usually derived by integrating the Newtonian 
stress-energy conservation equations over the interior of 
the body. Here, however, we do not assume the validity of 
the Newtonian equations in the interior. Equation 1)4. 2b|) 
contains all the information one needs in order to derive 
the explicit coupled equations of motion for the center of 
mass worldlines of each body in an A^-body system. We 
will describe this derivation later in Sec. IVII 

At post-l-Newtonian order, the formulae analogous to 
Eqs. (|4.2a|l and l)4.2b|l are 



S, = 



OO ^ 

^ [(; + 1) "Ml '-Gl + I "Ml °Gl 



1=0 



OO 

1=0 



"Ml "G.l + "Ml ^"G^l 



(4.3a) 



-SlHil + , , „ £tjfc "MjLHkL 



1 + 2 



l + l 



eijk "MjLHkL 



"([+'2)^ ^ ^i'' 3^ '^^ 

21^ + 7P + 151 + 6 



{l + l){2l + 3) 

OO ^ 

J2 n'^ik "^^^ 



1 + 2 

"M.L ^L - + - 



(^ + 1)' 



; + l 



kL- 



(4.3b) 
(4.3c) 



1=0 



For bodies in which the post-Newtonian field equations are valid everywhere, DSX derived the formulae (|4.3a|l - 
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H4.3c|l by using Newtonian^^ and post-Newtonian stress- 
energy conservation in the interior of the body In 



this section we will derive Eq. (|4.3b|l from the assump- 
tions listed above. The formulae (|4.3a|) and 14.3c|) giving 
the time evolution of the mass monopolc and the spin 
will be derived from the same assumptions in the second 
paper in this series W^. 

Note that the formulae (|4.3a|l - 14.3c|l are valid for 
all coordinate systems satisfying the assumptions listed 
above, not just for the body-adapted coordinate system 
discussed in Sec. IIII E)l above. 



B. Method of derivation 

We now turn to a description of the surface-integral 
method of derivation that we use. For this description 
we return, temporarily, to the context of the full, nonlin- 
ear equations of general relativity. The method is well- 
known and is described in Landau and Lifshitz |5lj| and 
in Misner, Thorne and Wheeler It has been pre- 

viously applied to the derivation of laws of motion by 
Thorne and Hartle |34l] . 

The method starts by fixing a coordinate system x'^ = 
{t,x^), and by writing Einstein's equations in that co- 
ordinate system in a form involving pseudotensors and 
partial derivatives as 



,a/3 ~ 



(4.4) 



Here g = det g^j,^ and the tensor density Tif^^^f^ is given 
by 



where 



a/3 



(4.5) 



(4.6) 



In Eq. 1)4.4(1 , T^'' is the stress-energy tensor and the pseu- 
dotensor T^" is given by 



167r 



,A0 ,Ai 



+ \ (2ff"^/^ - g^Pg^^) {2g, 



P0 ,A 



(4.7) 



V whose topology is that of the interior of a sphere, and 
whose boundary S = dV has spherical topology. We 
assume that the stress energy tensor T^" vanishes in a 
neighborhood of S. We define a quantity associated 
with S and with the choice of coordinate system by 



Phit) 



1 

16^ 



(4.8) 



Here d?^j is the natural surface element determined by 
the flat metric (dx^)^ -I- {dx'^)'^ + {dx^ f. Using the flat- 
space Gauss's theorem, the field equation 14.4|l . the sym- 
metry properties of Ti.'^"'^^ , and the assumption H4.1|l . the 
definition H4.8|l can also be written as 



PUt) 



d^x [(-g)r°' + r"] 



(4.9) 



Here d^x is the volume element associated with the flat 
metric (dx^)^ + {dx'^Y + {dx^Y . As is well known, for 
an asymptotically flat spacetime the quantity yields 
the ADM 3-momentum in the limit where S tends to 
infinity. For finite E, however, P|, does not have any 
invariant physical or geometric meaning. Nevertheless, 
we can still use this quantity as an intermediate tool in 
mathematical calculations in deriving relations between 
quantities whose transformation properties and meaning 
are well understood, such as the multipole and tidal mo- 
ments discussed in the previous section. For convenience, 
we will refer to P^ as the "enclosed 3-momentum" in E, 
even though it is not an invariant quantity. 

It follows from the formula H4.9|l and from the form 
[{~9)P^'^ + T^''] = of stress-energy conservation that 
the time derivative of the enclosed 3-momentum is 



pi 



r^d^Y., 



(4.10) 



The core of the surface integral method is to compute 
both sides of Eq. (|4.10(l explicitly. Namely, we compute 
the left-hand side by performing the surface integral 1)4. 8|) 
and taking a time derivative of the result, and we com- 
pute the right-hand side using Eq. (|4.7I) . These surface 
integrals are evaluated using the general solutions (|3.5a|) 
- H3.5c|) to the vacuum post-l-Newtonian field equations. 
The dependence on the surface S drops out and one ob- 
tains in this way relations between the moments and their 
time derivatives that reduce to Eqs. (|4.2b|l and l|4.3b|) at 
Newtonian and post-Newtonian order, respectively. 



This formulation of Einstein's equations is due to Landau 
and Lifshitz |5ll |. 

Next, in a surface of constant t within the domain of 
the coordinates, we consider a three dimensional region 



As is well known, the derivation of the spin evolution equation 
I4.ijcl requires only Newtonian-order stress-energy conservation. 



Note that the derivation of this conservation law does not require 
the assumption that Einstein's equations are valid for r < r_, 
since we can take the point of view that the stress-tensor T^"" 
is defined by Eq. 14.41 . In other words, a different theory of 
gravity could be applicable in the strong field region r < r_ , with 
the correction to the field equations being incorporated into the 
definition of T'^" . Our application of the conservation law I4.1UI 
to derive the equation of motion 14.3bt will therefore apply to any 
theory of gravity for which the vacuum field equations coincide 
with those of general relativity. 
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C. Post-2-Newtonian field equations 



and 



In applying the surface-integral method to derive the 
laws of motion to some post-Newtonian order to, one 
needs to compute some pieces of the metric to post- 
Newtonian order m + 1. This was emphasized in the 
original work of Einstein, Infeld and Hoffmann 0. For 
example, the Newtonian mass currents, whose conserva- 
tion law can be used to derive laws of motion at Newto- 
nian order, are the source for the gravitomagnetic poten- 
tial C [cf- Eq. H2.9c|l ]. Therefore, if one wants to avoid 
dealing with the matter distribution itself and to use in- 
stead only the far field metric to derive Newtonian laws 
of motion, one needs to use the post-l-Newtonian grav- 
itomagnetic potential. Note however that knowledge of 
the post-l-Newtonian scalar potential ip is not required 
for this purpose. 

Therefore, for our goal of computing the post-l- 
Newtonian laws of motion, we need to compute some 
pieces of the metric to post-2-Newtonian order. In this 
subsection we derive the harmonic gauge vacuum field 
equations satisfied by these post-2-Newtonian fields. 

It turns out that we need the post-2-Newtonian cor- 
rections only to gij and to ggi. We can therefore param- 
eterize the metric as the post-l-Newtonian metric (|2.4f) 
together with the appropriate correction terms: 

+ [2e^Q + 2e*(2$C» + + 0{e^)] dx'dt 

+ [S^j - 2e^^S^J + e*[(2$2 ^ + xkk)6^J - xy] 

+0{e'^)]dx'dx\ (4.11) 

Here the post-2-Newtonian fields are the vector field £,i 
and the spatial symmetric tensor Xij- We have chosen 
the particular parameterization of the post-2-Newtonian 
pieces of the metric (|4.11|l to simplify the gothic metric 
and the field equations. From the definition l|4.6(l of the 
gothic metric we obtain 

3'' = ^S'^ +e\,,+Oie^). (4.12) 

Inserting these expressions into the harmonic gauge con- 
dition ()2.5|l yields the equations 



(4.16) 



and 



r.,-Xfcfc+4V'-16$$ = 0. 



(4.13) 



(4.14) 



The harmonic-gauge vacuum field equations for the fields 
and x*-' can be derived by substituting the expansion 
(I4.12|l of the gothic metric into Eq. H4.4|l . The results are 



(4.15) 



We next compute the expansions of the enclosed mo- 
mentum and the spatial components T*-' of the 
Landau-Lifshitz pseudotensor that appear in conserva- 
tion law H4.1()|l . Substituting the gothic metric compo- 
nents (|4.12|) into the definition (|4.8|) of enclosed momen- 
tum yields 



where the Newtonian piece is 



1 



"PI = (t d,Cd'^, 



and the post-Newtonian piece is 



(4.17) 



(4.18) 



(4.19) 



The corresponding expansion of T*^ is, from Eqs. (|4.7|1 . 
fllHl and jUISll, 



(4.20) 



where the Newtonian piece is 

■T^^ = ^ (^d,^d,<i> ~ i<5,,afe$5fe$j (4.21) 
and the post-Newtonian piece is 



1 

47r 



d[iCk]d[kCj] + 2%$4) + 29(»*9j)V' 



4>2 



+2^k^^k^p + 3$ 



(4.22) 



D. Newtonian order derivation 

To illustrate the method of computation, we first derive 
the Newtonian law of motion Ij4.2b|) from the conserva- 
tion law <|4.10|l . A similar derivation has been given by 
Futamase 21J. We choose the 2-surface E to be the coor- 
dinate sphere r = R, for some R with r_ < i? < r_|_, and 
we henceforth drop the subscript S in P^ for simplicity. 
The formula H4.18|l then reduces to 



7-.2 r 

'P' - — $ n^djCiR-^) dn, 
IGtt J 



(4.23) 



v^r = C-i-i2<i>a,<f + 8%Cfc]5fc$ 



where the integral is over the unit sphere and — x-' /\:x.\. 
Plugging in the general solution (|3.8|) for the gravitomag- 
netic potential and using the identities ljA12p and (jA13|) 
yields 
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Zihdji 



1 ; 



.i=0 



L-l 



dr2 



-I r=R. 



E 

.i=0 



(-1) 



i+1 



°° f-lV+i 1 / 



L;=o 



(4.24) 



Since Z;^ and 1"-^ are STF on L, only the I — terms and the integrals ljA14p — (jA17p . First we write the 
can contribute [cf. Eqs. (|A14|I and HA15|l ]. This implies derivative of the expansion H3.5a|l of the Newtonian po- 



that 



tential as 



^ 16^ i?2 ^» 



(4.25) 



(4.27) 



1=0 



where we have used the relation (|3.1U|) . Note that the 
contribution from the gravitomagnetic tidal terms van- 
ishes identically. Thus, the Newtonian enclosed momen- 
tum is the time derivative of the Newtonian mass dipole. 

We turn next to the surface integral of the Landau- 
Lifshitz pseudotensor on the right hand side of Eq. I|4.10|l . 
To Newtonian order it is given by 



where 



V, 



* (2^ + 1)!! 



'Mr 



(4.28) 



and 



3 ' 

(4.26) 



(2^ + 1)!! „ 



r.i-l-3 



(4.29) 



from Eq. H4.21|l . This surface integral can be computed Substituting this into the right-hand side of Eq. H4.26|) 
using the expansion (|3.5a(l of the Newtonian potential and using the integrals ljA16|) and IjAlTp gives 
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i2l + 3)\\ 
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OO ^ 



(4.30) 



r 



Here the last line has been obtained using the expressions and r, and performing some easy algebra. We thus have 
(lOI^ and for and V^f in terms of "Ml, "Gl 
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the key result: 



•' 1=0 



(4.31) 



Substituting Eqs. 14.25|l and H4.31|l into the expansions 
(14.171) and (|4.2UI) and then into the conservation law 
(I4.1U|) finally gives the Newtonian law of motion l|4.2b|l . 



E. Post-Newtonian order derivation 

We now proceed with the derivation of the post-1- 
Newtonian law of motion (j4.3bp from the conservation 
law H4.10|l . Our method of derivation will be somewhat 
different from that used above in the Newtonian case. 
We start by describing the differences. 

The differences are related to the pieces of the com- 
putation that go to one higher post-Newtonian order 
(post-l-Newtonian order in the last subsection, and post- 
2-Newtonian order here). In the Newtonian case, we had 
available the explicit parameterization Ij3.5b|) - H3.5c|l 
of the post-l-Newtonian potentials in terms of (i) the 
Newtonian order moments "Ml, XJl; (ii) the post-l- 
Newtonian order moments ^"Ml, ^"Gl, Hl and Sl, and 
(iii) the gauge moments fiL and vl- The computation 
of the Newtonian enclosed momentum H4.18|l did involve 
the post-l-Newtonian potentials, but using the param- 
eterization H3.5c|l we found that the dependencies on 
the post-Newtonian multipole, tidal and gauge moments 
dropped out. Thus, we obtained equations of motion 
in terms of the purely Newtonian variables. Similarly, 
in the post-l-Newtonian case, the computation of the 
post-l-Newtonian momentum (|4.19ll involves the post- 2- 
Newtonian potentials and Xij- Those potentials can 
presumably be parameterized in terms of Newtonian and 
post-l-Newtonian moments, a set of post-2-Newtonian 
moments, and gauge degrees of freedom, via expansions 
analogous to Ij3.5bp - (|3.5c|l . Using such expansions we 
could in principle proceed as in the Newtonian computa- 
tion. However, this approach turns out to be extremely 
tedious^'', and we will proceed instead as follows. 

Our strategy will be to argue indirectly that there is 
no dependence on post-2-Newtonian degrees of freedom 
in the post-l-Newtonian momentum H4.19|l . The space of 
post-2-Newtonian potentials (^*, Xij) that satisfy in T> the 



The surface integrals encountered in the derivation of the evo- 
lution laws I4.ijal and I4.iicl for the mass monopole and spin 
are significantly simpler than those for the dipole evolution law 
14.3bt . We have derived explicit expressions for one particular 
solution of the post- 2-N ewtonian field equations for ^* and x'^ i 
and in a later paper l37l those expressions will be used to derive 
the mass monopole and spin evolution laws I4.3al and I4.3cl 
using the same explicit method as used here in the Newtonian 



field equations H4.15() and 14.16|l and the gauge conditions 
H4.13(l and 1)4.14(1 has the structure of an affine space. 
Given a particular solution (CpiXp )j ^my other solution 
can be written as the sum 



f = Cp + CJ; 



= Xp'+Xh' 



(4.32a) 
(4.32b) 



of the particular solution and a homogeneous solution 
(^h' Xh )' where the homogeneous solution satisfies homo- 
geneous versions of the field equations H4.15|l and ((4.16|l 
and gauge conditions (|4.13ll and 14.14|1 : 

= 0, \/\^ = 0, (4.33a) 
5,Xh^ = 0, ~ Xf = 0. (4.33b) 

The general solutions to Eqs. H4.33a|l can be expanded as 
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X,i9i—--VK,LX^, (4.34a) 



1=0 



1 1 



n ^^^^^^ W " i!^'^^^^ ■(4-34b) 



These multipole expansions define the multipole and 
tidal moments Xi^, WiL, C'ijL and Bij^, all of which 
are STF on the indices L. Inserting these expansions 
into the gauge conditions 14.33b|l yields the constraints 



Ci<L> — 0, 
IX^L> + CjjL = 



B^JJL = 0, (4.35a) 



Next, we insert the decompositions (|4.32a|l - (|4.32b|) 
and the expansions ((4.34a|l - ((4.34b|) into the formula 
(|4.19(l for the post-l-Newtonian momentum. Using the 
integrals HA14|) ~ (|A17p then gives 



1 1 ■ ■ 



(4.36) 



Here as before we have chosen the 2-surface E to be the 
sphere \x\ — R. The last three terms in this expression 
give the dependence of ^"P* on the homogeneous solution. 
From the constraints ((4.35a(l - ((4.35b(l it follows that the 
sum of these three terms vanishes, so that 



1 

16^ 



f [d,(l + x;']d'^,. (4.37) 



Thus, the enclosed momentum '"'Pi is independent of 
which solution of the post-2-Newtonian field equations 
is chosen. Therefore, it must be a function only of the 
Newtonian and post-l-Newtonian fields, or equivalently 
of the Newtonian and post-l-Newtonian multipole and 
tidal moments, as well as of the radius R. 



27 



Next, we show that the enclosed momentum can be 
written as the sum of the time derivative of the post- 
Newtonian mass dipole p"Mi and terms that are inde- 
pendent of '■"Mi. To see this, note that any solution 
(<&,(■%''/') of the post-l-Newtonian field equations (|2.9a() 
- (|2.9c|l and gauge condition (|2.6() can be decomposed as 



$ = $0, 



C 
'"MA 



Co, 

1 



(4.38) 
(4.39) 



Here ($o, Co, V'o) is another solution with vanishing post- 
Newtonian mass dipole. Inserting this decomposition 
into the post-2-Newtonian field equations and gauge con- 
ditions (|4.13|) - H4.16|l yields a corresponding decomposi- 
tion of the post-2-Newtonian potentials: 



x'' - xo^ 

4 - 4o 4 M 



(4.40) 
(4.41) 



Here the potentials (<&o, Co, V'o, Xo , ^o) ^'^'^ ^ solution of 
the field equations and gauge conditions and arc indepen- 
dent of '■"Mi. Inserting the decompositions (|4.40() and 
(|4.41(l into Eq. H4.19|l gives 



1 
16^ 



(f^j. (4.42) 



Inserting this into the expansions H4.17|l and H4.20|l and 
then into the conservation law (|4.10l) now yields 



(4.43) 



Equation (|4.43|) is the key result of this subsection. 
It shows that the second time derivative of the post- 
Newtonian mass dipole can be expressed purely in terms 
of the Newtonian and post-l-Newtonian fields, and is in- 
dependent of the post-2-Newtonian degrees of freedom. 
This independence was derived above for the first two 
terms on the right hand side of Eq. H4.43|l . and for the 
third term it follows from the fact that the expression 
(|4.22l) for p'T*-' depends only on the Newtonian and the 
post-l-Newtonian fields. Therefore the right hand side 
of Eq. (|4.43|l is some function of the surface S as well 
as of the multipole and tidal moments of the solution 
(<&, CNV') [since the multipole and tidal moments of the 
solution ($0, Co, V'o) coincide with those of the original 
solution {^,(\'>P) except for the post-l-Newtonian mass 
dipole] . 

We can deduce some properties of the functional de- 
pendence of p°Mi on the moments as follows. From the 
post-2-Newtonian field equations and gauge conditions 
()4.13(l - H4.16|l it follows that the potentials and Xo 
depend linearly on Q and "00, and quadratically on $0, 
$0, and 5[*Co'- Also from Eq. it follows that ^'T*-'' 



depends quadratically on $, $, d^'^(^\ and C* + d''ip, all 
of which are independent of the gauge moments /i^ and 
h'L- From Eq. 14.43|l and the expansions (|3.5a|l - H3.5c|l 
it now follows that 



TA --Ml, "Ml. "Ml. 'XJl, "Gl, "Gl, Hl, Hl, 



Sl.Sl, '"Ml, ''XJl'tR 



"Ml, "Gl, Hl, Sl, 'f^L, v'l; R 



(4.44) 

Here Qi is a linear function of all of the moments that 
appear as its arguments, and can be an arbitrary function 
of the radius R that defines the 2-surface S. Similarly 
the function Ti is a quadratic function of all the moments 
that appear as its arguments, and can be an arbitrary 
fimction of i? 

In Appendix lUl we compute explicitly the linear term 
and show that it vanishes identically:^® 



e» = 0. 



(4.45) 



Next, the left hand side of Eq. (|4.44|l is independent of the 
radius i?, as are the definitions of all the moments that 
appear as the arguments of the function J^i. It follows 
that Ti is independent of R, as one would expect. Using 
these simplifications we can rewrite Eq. (|4.44ll as 



TA "Ml, "Ml, "Ml, "Gl, "Gl, "Gl,Hl,Hl, 



(4.46) 



Sl,Sl, ""Ml, '"Gl 



where the function J-i is now a quadratic function of all of 
its arguments. This quadratic function could in principle 
be computed from the expression H4.43|l . However, it 
is simpler to appeal to a special case from which the 
functional form of J^i can be deduced, as suggested in a 
different context by Thorne and Hartle [s^l • 

Specifically, we now specialize to the case considered 
by DSX where the post-l-Newtonian field equations are 
assumed to hold throughout r < r_ . Our analysis applies 
to that special case, and therefore the functional J^i coin- 
cides with that computed by DSX, given in Eq. (4.21b) 



It is easy to see that the right hand side of Eq. i4.44i must be 
independent of the radii r_ and r+ that define the domain T). 
The faet that the right hand side of Eq. 14.441 is independent 
of the gauge moments fi^ and ul can alternatively be derived 
as follows. As discussed in Sec. 1111 Dl above, by making a gauge 
transformation of the type I2.22i we can alter the values of the 
gauge moments fiL and ul without altering any of the tidal and 
multipole moments "M^, =Gi, ^"Ml, ^"Gl, Sl, Hl or their 
time derivatives. Under such a transformation, the left hand side 
of Eq. I4.44i is invariant. It follows that the right hand side does 
not depend on the gauge moments or their time derivatives. 
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of Ref. [lal and in Eq. Ij4.3b|l above. The derivation of 
the form of Ti in this case is reviewed in Appendix ^ 
The key point here is that the general argument of this 
subsection estabhshes the result (I4.3b|) up to the values of 
the coefhcients of the terms on the right hand side, and 
our argument shows that those coefficients are univer- 
sal, applying both to the case of weakly self-gravitating 
bodies analyzed by DSX and to the case of strongly self- 
gravitating bodies considered here. 

This argument, which enables us to avoid doing the 
surface integral H4.43|l explicitly, could of course also be 
used to avoid doing the surface integrals in Eq. (|4.2()|) 
when computing the Newtonian laws of motion. In that 
case, however, we were able to check explicitly that the 
surface integral method gives the correct answer. 

Similarly, we could deduce that the linear term Qi 
must vanish by comparison with the case of weakly self- 
gravitating bodies. Therefore the explicit verification of 
this result in Appendix lUl is not really necessary. That 
verification is useful, however, as a consistency check of 
our argument and formalism. 

Thus, we have established that the law of motion (|4.3b|) 
is valid not just for the class of weakly self-gravitating 
bodies considered by DSX, but also for the more general 
class of strongly self-gravitating bodies considered here, 
subject to the assumptions outlined in Sec. II V Al 



V. AN N-BODY SYSTEM: FOUNDATIONS 

We now turn to an analysis of a system consisting of 
N bodies with arbitrary internal structure. The bodies' 
masses can be comparable, but their typical separations 
must be large compared to their masses in order that 
their gravitational interactions are well described by the 
post-l-Newtonian approximation^^. In this section we 
lay the foundations for our analysis by defining local coor- 
dinate systems associated with each body, and an overall 
global coordinate system. We also derive relations be- 
tween the moments that characterize the potentials in 
each of these coordinate systems. In Sec. I VII below we 
will combine the results derived here with the single-body 
equation of motion Ij4.3b|l derived in the previous section 
to obtain the explicit form of the Af-body equation of 
motion. 



A. Assumptions 

We start by describing our assumptions. Wc con- 
sider a system of N bodies, labeled by the index A with 



\ < A < N . We associate with each body a world tube 
Wa containing the region where the stress-energy tensor 
is nonzero, and also containing the strong-field region as- 
sociated with the body. Our key assumption is that the 
post-Newtonian equations are satisfied everywhere out- 
side all of the worldtubes YVa- 

More precisely, we make the following assumptions: (i) 
For each A there exists a coordinate system {sat'u'a) 
the type discussed in Sec. IIIIDI above which covers the 
Ath body. Thus, there exist radii t-^a and r+,A such 
that the range of the coordinates includes the product of 
the ball \yA\ < t+,a with an open interval of time, and 
that the coordinates are harmonic, conformally Cartesian 
and body-adapted in the buffer region 



Ba = \{sa,V^j^ r-^A<\yA\<r+,A} ■ 



(5.1) 



(ii) The various buffer regions Ba are non-intersecting 
(see Fig. ^ above). (iii) We define the world 
tube associated with the Ath body to be = 

|(sA,y;ji) Iz/aI < and we define the spacetime re- 

gion T) to be the complement of the union of all the world- 
tubes, 

N 



I? = X \ y wa, 



A=\ 

where M. is the entire manifold. The vacuum Einstein 
equations are satisfied for the one-parameter family of 
metrics H4.11|l on the spacetime region T). (iv) There 
exists a conformally Cartesian and harmonic coordinate 
system {t,x^) which covers all of V. We will call the 
coordinate system {t,x^) the global frame, even though 
it does not cover the entire manifold. We will call the 
coordinate systems {satVa) ^^'^ body frames. 

Note that in the context of the spatially non-compact 
domain T), the meaning of the 0(e") symbols that 
appear in Eqs. (|2.2() and H4.11|l corresponds to pointwise 
convergence, not uniform convergence. As is well known, 
solutions of the post-l-Newtonian field equations are 
not good approximations to exact solutions at distances 
> 1/e. That is, although they work well in the near zone 
they break down in the local wave zone 38] . In order to 
obtain solutions which are good approximations every- 
where one has to perform a matching of post-Newtonian 
solutions onto radiation zone post-Minkowskian solu- 
tions; see for example Blanchet 3j and references therein. 
However, the corresponding corrections to the near zone 
gravitational fields and to the dynamics of the bodies 
arises at post-2. 5-Newtonian order, and can thus be 
neglected for the post-l-Newtonian analysis of this paper. 



Our formalism and derivation does not require that the bodies' 
separations be large compared to their typical sizes. However, 
that requirement is in practice necessary if one wants to achieve 
good accuracy using a truncated version of the equation of mo- 
tion 16. Hi containing only a small number of multipoles. 



B. Body- frame multipole and tidal moments 

In each local coordinate system (s^i, u^'a) we define mul- 
tipole and tidal moments "Aff(sA), '°g£{sa), ^''M^{sa), 
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^"Gf{sA), S^{sa) and H^{sa) according to the prescrip- 
tion described in Sec. 1111 iH above. We have added super- 



scripts A to these moments to denote the ^th body. The 
corresponding expansions of the potentials are 



/=o ' \yA\ ■ 



1=0 



^ ' "Mf{sA)dL- 



vaI n 

"Gti.sA)y^}. 



CO 

(Hsa.v'a) = E 



1=0 
1 



qA 



dr 



\yA\ 



Y^-^^ji<aiH^_l^j{sA) 



4(2? - 1) ^A , ^e 
2^ ^ 2^ ^<L-l\SA)Oai>t 



1=0 



IvaI IV 
r 



(5.2a) 



(5.2b) 



(5.2c) 
(5. 2d) 



where overdots mean derivatives with respect to the time 
argument. These expansions are obtained from the ex- 
pansions H3.28all — H3.28c(l by replacing with and 
t with SA, by adding superscripts A to the potentials 
and the various moments to denote the Ath body, and 
by omitting the gauge moments /il and which vanish 
since we have specialized to body-adapted gauge. The 
specialization to body-adapted gauge also implies that 



0, 



(5.3) 



cf. Sec, rrrrnl above. 



C. Configuration variables for tlie Atli body. 

For each body, there is a non-empty region of over- 
lap between the domain of the body-frame coordinates 
{satV''a) ^^'^ domain of the global- frame coordinates 
(t, x*), namely the buffer region Ba defined by Eq. H5.1|l . 
Both of these coordinate systems are harmonic and con- 
formally Cartesian, and therefore the mapping between 
the two coordinate systems can be parameterized using 
the general analysis of Sec. Ill Bl above. From Eqs. (|2.17|) 
- (|2.18c|l it follows that there exist functions z/^(sa), 
a^(s^), hf^{sA) and R^{sa) and a harmonic function 
/3f^(s^,y^) such that in Ba 



y\ + zf{sA 



)+e'hf{sA,y'A)+0{s'), 



SA 



'a^{sA,y\)+e^f3^{sA,y\) + 0{e'^), 



(5.4) 



where 

= af{sA)+y{zf{sA), 
hf = h^^isA)+e,jky'ARt{sA) + -zfisA)yWA 



(5.5) 



1 



-y\d^{sA) - yWAifisA) + -y\zf{sA)if{sA) 
1 

"2' 



'-zf{sA)zf{sA)y'A 



(5.6) 



and 



vWa 



^'^'i{sA)y\ + \a^{sA) 



-f3^{-SA,y\). 

(5.7) 



Because the body-frame coordinates are uniquely de- 
termined, the functions zf{sA), a^(syi), h^^{sA) and 
RkisA) acquire the role of configuration variables that 
specify the location, orientation etc. in the global coordi- 
nates (i, a;*) of the local rest frame attached to body A. 
This role is in contrast to the role of the corresponding 
variables in Sec. above, which were freely specifiable 
functions. The task of determining the motion of the N 
different bodies reduces to solving for the time evolution 
of these configuration variables. Below we will show that 
some of these variables can be obtained by solving dif- 
ferential equations, and the remainder are obtained from 
algebraic relations. 

The particular combination of these configuration vari- 
ables that enters into the equation of motion which we de- 
rive below is the center of mass worldline. In the special 
case where the post-l-Newtonian equations are assumed 
to hold inside each body, this center of mass worldline is 
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defined simply as the origin of spatial coordinates of the 
body-adapted^^ coordinate system (satVa)- "^^^s world- 
line can be expressed in terms of the global-frame coor- 
dinates in parametric form as 

x'{sa) = zf{sA)+s^h^^{sA) + 0{e*), (5.8) 
t{sA) = SA+e^a^isA) + Oie^), (5.9) 

from the coordinate transformation H5.4(l - H5.7|l . Elimi- 
nating sa gives = "'^zf{t), where 

-zfit) ^ zfit) + [h^S) " if(.t)atit)] + 0{e% 

(5.10) 

Here the superscript "cm" means "center of mass" . 

Consider now the more general context where the post- 
1-Newtonian equations are not assumed to hold inside 
each body. Then, the body-frame coordinates {satTJ^) 
can be arbitrary for \yA\ < so the worldline in 

spacetime of the origin ua = of these coordinates has 
no special significance. Nevertheless, we can still use Eqs. 

' to define the function •"^zf{t). That IS, we 

define "''^zf-{t) to be the image of the origin = not 
under the true coordinate transformation, but under the 
extension to \yA\ < T-^a of the formulae (|5.4|) - (|5.7() 

More precisely, of the slightly modified body-adapted coordinate 
system discussed in the last paragraph of Sec. lIII Cl above. whose 
domain of definition includes the interior of the body. The dif- 
ference between this coordinate system and the body-adapted 
coordinate system arises only at order 0(£*) in the time coor- 
dinate, which does not affect the definition of center of mass 



tp^{t, X^) 



where 



Here the superscript "g" on the potentials and on the 
moments stands for " global" . 

The form of the expansions H5.11a|l - (|5.11c|l is dictated 
by the following considerations: (i) The expansions take 
the form of a linear superposition of solutions, one for 
each body A. This follows from the linearity of the vac- 
uum field equations (|2.9a|) - (|2.9c|l . (ii) We choose to 



which a priori are only valid for r-^A < \yA\ < f+,A- 
The resulting function ^"'"zf{t) continues to characterize 
the location of the local rest frame attached to body A, 
even though it no longer corresponds to a worldline in 
spacetime^^, and even though the location = °'^zf{t) 
will in general be outside the domain of definition of the 
global coordinates. We will continue to call this function 
the center of mass worldline, in a slight but conventional 
abuse of terminology. 



D. Global- frame multipole moments 

In this section we define, for each body A, multipole 
moments associated with the global coordinate system 
(i, a;'). We define the global- frame multipole moments 
"Af|'^(i), ^-Mf^it), Sl^{t) and nl^(t) to be the mo- 
ments about the Newtonian-order center-of-mass world- 
line X = (t) of body A [cf . Eq. (|5.10l) above] , using the 
prescription discussed in Sec. lHlEl Using these multipole 
moments we can write down multipole expansions of the 
global-frame potentials, which we denote by ($s, (^s, ^/jS), 
that are valid on the entire domain V: 

worldline to post-l-Ncwtonian accuracy. 

The function ''"'zf-{t) does however transform like a worldline 
under the group 12.171 of post-l-Newtonian coordinate transfor- 
mations. 



I 

use a gauge for the global coordinates in which all the 
potentials go to zero as |a;| — > oo. This eliminates any 
tidal terms associated with acceleration of the reference 
frame, cf. the discussion in Sec. IH Bl above, (iii) There 
are no other tidal terms, since the terms in the sum over 
B with B A play the role of tidal terms for body A. 
With this identification, the expansions (|5.11a|l - H5.11c|) 
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agree with the formulae H3.58a|) - (|3.58c|l of Sec. IIIIEI 
above that define the multipole and gauge moments, in 
the buffer region Ba about the Ath body. 

We will derive transformation laws relating the global- 
frame multipole moments "M|'^, p°M|'^ and 5*1'^ to the 
body- frame multipole moments "M^, '"'Mf and Sf in 
the next subsection. Note that the moments that would 
be measured by observers residing in the buffer region 
Ba about the ^th body are the body-frame moments 
and not the global-frame moments. 

We next discuss the gauge freedom in the global coor- 
dinate system. If we make a gauge transformation of the 
form 12.22|l with the harmonic function /?h chosen to be 
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(5.12) 



then the gauge moments ii^^{t) transform according to and 



s,A _ „g,A ^ a + l)(2/ + 3) 
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cf. Eq. (|3.54f) above. Therefore there is enough freedom 
to set 



Mi'^ = 



(5.14) 



for all A and for all / > 0. This requirement, together 
with the requirement that the potentials go to zero as 
I a; I oo, reduces the residual gauge freedom to the post- 
Galilean transformation group discussed in Sec. IIIBI 



E. Computation of body-frame tidal moments 

Our goal is to deduce equations of motion for the A^- 
body system from the single-body equation of motion 
(|4.3b|) . To this end, we would like to compute the body- 
frame tidal moments "Gf , ^'Cf and felt by body 
A in terms of the body- frame multipole moments "Mf , 
p°M£ and and also the configuration variables of the 
other bodies B with B A. We shall perform this com- 
putation in stages, by relating both sets of quantities to 
the global-frame moments. 

We start by expanding the time derivatives that appear 
in the expansion IjS.llbp of the global-frame potential 



and by expressing the results in terms of STF tensors. 
Using the gauge specialization 1)5.14(1 this computation 



gives 
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where the STF tensors Nfj^ and Pl'^ are given by [cf. 
Eqs. (|3.61b|) and I3.61cl) above] 
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Next, we expand the global potentials in the buffer 
region Ba of body A using the Taylor series 
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T^{z)^id^\z~xn^^,. 
For p = —1 we have 
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Substituting Eq. (|5.18() into Eqs. H5.11a|l . H5.11c|l and 
((5.15(1 and using the identity 
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yields 
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where 
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compute the body-frame moments in terms of the global- 
frame moments, as in Sec. IIII CI above. At Newtonian 
order we obtain [cf. Eqs. - (|3.35b|l above] 
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=Gf = TJl'^-nAf, (5.25b) 



where the nonzero inertial moments A* are given by 

1 



and 
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(5.26a) 



(5.26b) 



Transforming the gravitomagnetic potential gives [cf. 
Eqs. - (|3.36b|) above] 



(5.27a) 



Here "G*!'^, y^^"^, and Jf^^ are global- frame tidal 

moments. The post-Newtonian moments Ff^'^ , ^uf^ ^-'^'i 
Jf^^ could be parameterized in terms of the irreducible 
global-frame tidal moments ""Cl'"^, Hf^^ , and if de- 
sired via equations analogous to Eqs. 13.61dp . H3.61e|l 
and H3.58e|l . Here, however, it will be more convenient 
to work directly with the moments Ff^' , and J|' . 
The function z^{sb) that appear on the right hand sides 
of Eqs. (|5.23al) - H5.23dp is evaluated at = <. 

Next, we apply the coordinate transformation 15. 4|) to 
the global potentials (|5.22a|l — (|5.22c|l using the formulae 
(|2.2Ua|l - H2.20c|) . We parameterize the harmonic func- 
tion f3^ which appears in Eq. H5.7II as 

°° I'-lV+i 1 



1=0 



00 



(5.24) 



where the tensors and are STF. Comparing the 
result with the expansions H5.2a|l - (|5.2c|l allows us to 



y^^ = - 4i,^ "G«/ - - /!Af^, (5.27b) 
where the nonzero inertial moments A^^ are given by 
A^ = h^^ + zfzfzf - t^jkzfR^ - 2a^z^, (5.28a) 



and 
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Next, by combining the transformation laws (|5.27a|l and 
(|5.27b(l with the decompositions of Z;^ and Y^^f^ given by 
Eqs. (|5.2c|) - H5.2dp . the decomposition IjS.lldp of Z: 
and the gauge condition (|5.14() . we can solve for the 
ordinate transformation functions Xf^ for I > and 



L ' 
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for / > 1. The result is 



(5.29a) 
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and ing the definition H5.16|l of 7V|' , the formulae (|5.29a|l 

and Ij5.29b|l for and t^, and the identities HA8|I - 
r^L = ~ 4ii^ "G|'^ - IIA%^^^. (5.29b) ^ gives 

Finally, matching the post-Newtonian potentials and us- 
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Here then nonzero inertial moments A^'' are given by 



^3 "c J 



Af- = 



~A uA 



-od'zt 



+Q!;^ifif, (5.31a) 



+ -i^z^z^ 



2 « J J 



(5.31b) and 
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In Eq. (|5.30bl) it is understood that the moments Y^'^ 

and a£ are zero for I — 0. 

The le ft-hand sides of Eq s. ^5.1^5a|l . ^5.25bp . (|?:77ajl . 
(|5.27b|) . (|5.30a|l and l)5.30b|l are functions of the time 
coordinate sa of the body-adapted coordinate system for 
the Ath body, cf. Eqs. H5.2a|l — H5.2c|) above. The right- 
hand sides are expressed as functions of sa by evaluating 
the global moments, which are functions of the global 
time coordinate t, at t = sa, cf. the discussion in the 
last paragraph of Sec. IIIICl above. 

Finally, by combining the transformation laws Ij5.25b|l . 
(|5.27b|l and (|5.30bp with the gauge speciahzations (|5.3(l of 
the body-adapted coordinates we can deduce the values 
of some of the configuration variables of the Ath body. 
We obtain 



= / dsA 
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(5.34) 



The only remaining configuration variables that are un- 
determined are the variables h^^ and zf that determine 
the center of mass worldline (|5.10() . 

To summarize, the main results of this subsection are 
the explicit expressions (|5.25bp . (|5.27b|l and (|5.30b(l for 
the body- frame tidal moments "Gf^, p"G^ and which 
act on body A in terms of the configuration variables of 
all the bodies, as well as the global-frame mass and cur- 
rent moments "M^' , ^"Mf^' and 5^' of the other bod- 
ies. These expressions are given by combining (j5.25b|l . 
(|5.27b|) and (|5.30b|) with Eqs. ()5.ij3a|l — (|5.23d|) . Also, 
the global- frame multipole moments "M|'^, ^"Mfj^ and 
5|' of body B can be reexpressed in terms of the body- 
frame multipole moments "A/^, ^"Af^ and of that 
body using the relations H5.25a|l . (|5.27a|l and (|5.30a|) . 
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F. Definition of body frame multipole moments 

Mt and St 



As discussed in the introduction, it is useful to use in- 
stead of the body-frame muhipole moments Mf and Sf 
a modified set of body-frame moments defined as follows. 
We define for each body A a coordinate system {s^,yf) 
which is identical to the body-frame coordinate system 
{s^, yf) except that it is non-rotating with respect to the 
global frame coordinates {t,x^) (i.e., non-rotating with 
respect to fixed stars). We define the moments Aif^{t) 
and Sf{t) to be the multipole moments of body A in 
this non-rotating coordinate system, expressed as func- 
tions of the global time coordinate t. These are given by 
the equations 



A a', I 



(5.35) 



and 



where SA(t) is the value of the body- frame time coordi- 
nate sa evaluated at what would be the intersection of 
the worldline of body A with the spacclikc hypcrsurfacc 
of constant t. From Eq. 15.4|l this function is given by 

SA{t)^t-e^a^{t) + 0{e'^). (5.37) 

Also the rotation matrices U^"" are defined by the for- 
mula 



U, 



A a' 



A J- e-2<r 



From Eqs. (E2I), (gTI) and 

these moments as 



(5.38) 
we can write 
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+ l<^jk<ai "Ml 

St + 0{e^). 



+ 0(£^). (5.39) 
(5.40) 



As indicated by Eq. (|5.4U|I . when working to post-1- 
Newtonian order we can identify the moments iS^ and 
Sf. Nevertheless it might be useful in some circum- 
stances to use the more accurate relation (|5.36|) . for ex- 
ample for systems which evolve for sufficiently long times 
that the rotation matrices ° become significantly dif- 
ferent from unity. 

All the tools are now set up to compute explicit 
equations of motion for the center of mass worldlines. 



VI. EXPLICIT EQUATIONS OF MOTION FOR 
AN N-BODY SYSTEM 



In this section we derive explicit equations of motion 
for the center of mass worldlines '""z/^(t) of each body 



as seen from the global coordinate system, by combining 
the single-body equations of motion (|4.2b|l and Ij4.3bp 
with the moment transformation formula derived in Sec. 
IVI above. 

We start by deriving the well-known Newtonian equa- 
tions of motion, in order to illustrate the computational 
method. The Newtonian single-body equation of motion 
(|4.2b|) applied to body A implies that 
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(sA, 



-GfLi^A) 



(6.1) 



since the body-adapted coordinates are mass-centered, 
i.e. "Mf = for all A. Using the relation H5.25bp be- 
tween the body frame tidal moments "Gf^ and the global 



frame tidal moments °G|'^ 



"Gf' 



we can rewrite this as 
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(6.2) 



Here the acceleration zf of the Newtonian-order center of 
mass worldline has appeared via the transformation law 
for -Gf. Next, we substitute the expression H5.23a|l for 
the global-frame tidal moments "Gf^* in terms of mass 
multipole moments °M|'^ of the other bodies, and use 
Eq. 15.25a|) . This gives 



EEE 

B^A k=0 1=0 



AY 



-Mi 



kill "M^ 



M 



K 



(6.3) 



"M^ are functions of sa, 



Here the quantities and 

while the quantities and "M^ are functions of ss, 
evaluated at sb — sa- Writing the dependent variable 
as t instead of sa and using the definition (|5.19() of 7] 



K 



we can rewrite Eq. H6.3I) in the more explicit, well-known 
form 



vf:f:i#^"Mi(, 
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kill "M^ 



(6.4) 



The analogous computation carried to post-1- 
Newtonian order is similar but much more involved. We 
start by focusing on the first two terms on the right hand 
side of the single body equation of motion (|4.3b(l . and 
evaluating explicitly the I = pieces using Eqs. (|5.25b|) 
and Ij5.30b|l . [As before the / = 1 pieces vanish since the 
body-frame coordinates are mass-centered, by Eq. H5.3|) ]. 
Using the Newtonian equation of motion H6.3|l . the def- 
inition H5.39|l of A^f and the definition (|5.1U|) of """zf, 
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the result can be written in the form 



Here by analogy with Eq. 



1=2 



we have defined 



jk<ai >-JL-l>j^fe 



(iii) Compute the global- frame tidal moments "G^^, 
Y^'^, Jf^^ and Ff^'^ of body A in terms of the global- 
frame multipole moments of body B using Eqs. (|5.23a() 
- (|5.23dp . (iv) Compute the body- frame tidal moments 

(6.5) , =Gf and "'Gf of body A in terms of its global frame 

tidal moments °C?|"^, Y^^Jf;^ and Here the re- 

sults are given by Eq. (j5.25bp for °Gf , and by Eq. H5.30b|) 
for ^Xj^; note that for the required values of Z (Z > 3) 
the last three terms in Eq. Ij5.30b|l do not contribute. For 
we have Hf = by the gauge condition (|5.3|l . while 

(6-6) for Z > 2 we obtain from Eq. H5.27b|l that 



and gf is defined to be all of the terms on the right 
hand side of Eq. (|4.3b|l except for the first two terms 
(with superscripts A added to all the moments) . Also we 
define 
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In order to explicitly evaluate the tidal moments that 
appear on the right hand side of Eq. I|6.5|) , we perform the 
following sequence of moment transformations: (i) Start 
with the body frame multipole moments "M^, ''"M^ and 
Sf^ of body B. (ii) Compute from these the global frame 
multipole mo ments ° Af s-^, N f; ^ Pf '^ a nd Zf ^ of body 
B usin g Eqs . (jrnajl . ^5.11d|l . (|3.61b|) . (|XHTc| . (jF^TZajl . 
(|5.29a|l and H5.30a|) . The results are 
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(v) By combining the preceding steps, all the moments 
can be expressed in terms of the multipole moments 
"Mf , p-Aff and S'f of body B. 



In the resulting expression, we eliminate the variables 
in favor of ""zf for all C using the definition H5.10|l . 
and we eliminate the moments "M^ , in favor of A4£, 
Sf for all C using the definitions and (|On|) . These 

substitutions generate correction terms only in the 0(6*^), 
Newtonian terms in Eq. (|6.5|) . and not in the 0{s^), 
post-Newtonian terms^°, since we drop all terms of or- 
der 0{e^). We also eliminate using Eq. H5.32|l . The 
resulting expression then depends only on the variables 
""zf , A42 and iS£; the dependencies on the variables 
i?^, and h^.^ cancel out. Lastly, we set to one the 
formal expansion parameter e. The result of this tedious 
computation is 



B^A B^AC^A 

+ E T.^t''''it) + 0{e% (6.11) 

B^A C^B 



where 



36 



In particular since the current moments Sf^ do not enter at New- 
tonian order, there are no correction terms generated when one 
eliminates in favor of S^. Thus, we are free to use either 5^ 



or in the equations of motion, to post-l-Newtonian order, as 
noted in Sec. IVn above. 
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(i)p^sc _ (-l)^+^+^/(2fc + 2; + l)!!(2p + 1)!! a^^^b^.^c 

^^^^ ~ fc!/!p!(2p+l)A^^ MlMj^Mp, (6.14a) 

(2) 25ABC _ (-l)Wi(^ + 4)(2fc + 2/ + l)!!(2p + 2g+l)!! ,.a,.b^Q ,.c 

- fcmp!g!(2fc + 2/ + l)M^ ^^^^MA^p, (6.14b) 

(3) 2,ABC _ (-l)Wi(^ + i)(2fc + 2; + l)!!(2p + 2g+l)!! . Al^^ 

(4)pABC _ (-l)''+^+^/(2fc + 2; + l)!!(2p + 2g+l)!! ^ . ^ ^4^) 

^iKLPQ - k\l\p\q\M^ M,lMj,^^Mp, (6.14d) 

(5)pABC _ 3(-l)'-+n2fc + 2^ + 3)!!(2p + 2g+l)!! A^^^^^ :K 

(D^^ABC _ (-l)'-+^+^(fc + l)(2fc + 2/ + l)!!(2p + 1)!! b^.^c .515 . 

^^^^ " kmp\{2p+l)M^ ^lM^Mp, (b.lSa) 



(6.15b) 



(3)pABC 4(-l)^+^(2fc + 2^ + l)!!(2p + 2q + 1)!! . 

^^^^•^ - fc!;!p!g!(2fc + 2i + l)M^ Mr.M^ ^^Mp, (6.15c) 
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(4)jyABC 



(_l)fc+P+i(2fc + 21 + l)!!(2p + 2q + 1)!! 



k\l\p\q\M^ 
(2fc+ l)(/c + 2 



2fc + 3 



1 



2fc + 2r 



>3 



(6.15d) 



iKLPQ 



{-lf+P{l + 2){2l + l)(2fc + 2^ + l)!!(2p + 2g + 1) 
kW.p\q\{2l + 3) 



(6.15e) 



^ijKLPQ 



where 



3(-l) Wi(2A: + 2; + 3)!!(2p +2q+ 1)!! ^ . >(g ^ 



fc!Z!p!g!7W^ 



M 



(6.15f) 



(6.16) 



Here we have denoted by the unit vector pointing 
from the center of mass worldhne of body A to that of 
body B: 



where 



rBA (t) 



rBA{t) 



(6.17) 



(6.18) 



We also have defined = and vf"^ 



12cll depend 



The right hand sides of Eqs. H6.11|l - 
on the time derivatives Sf, and of the bodies' 
spins and mass monopole moments. These dependencies 
can be ehminated using the single-body laws of motion 
(I4.2a|) . (|4.3a|) . and (|4.3c|) . the last two of which are de- 
rived for strongly self-gravitating bodies in paper II [33 . 
For the case of the mass monopole moments, this proce- 
dure generates terms that are of post-2-Newtonian order 
which can be neglected. Thus, all time derivatives of 
mass monopoles appearing in Eqs. (|6.11|) - (|6.12c|l can 
be neglected. In other words, we can make the following 
substitutions in Eqs. H6.13a|l - (|6.15f|l : 



Mi 



Mi 



{l-Si„)Mi. 



(6.19) 



For the case of the spin time derivative terms, using Eq. 
(|4.3c|l together with Eqs. (|5.23a|) and H5.25b|l we obtain 
modified values of the coefficients (|6.13a() (j6.15f|) which 
are listed in Appendix IFl 

A simple special case of the above equations of motion 
is the non-spinning point particle model, or monopole- 
truncated model. This is obtained by setting to zero all 
the mass multipole moments for / > 1, and all the 
current multipoles Sf. In this case Eq. (|6.11|l reduces to 
the well-known Lorentz-Droste-Einstein-Infeld-Hoffmann 
equations of motion , which were also reported in 
the first DSX paper [Eq. (7.20b) of Ref. 0]. 

A second special case is the spinning point particle 
model or monopole-spin truncated model, obtained by 



setting to zero all the mass multipoles "Mf and ''"Mf 
for Z > 1, all the current multipoles S*^ for / > 2, but 
allowing non-zero spins Sf. For this case our general 
equation of motion H6.11|l reduces to the equations of 
motion obtained for this case by DSX [Eqs. (6.30) — 
(6.34) of Ref. ^]. 

Finally, we can obtain an explicit expression for the 
angular velocity (|5.34l) parameterizing the dragging of 
inertial frames by using the Newtonian equation of mo- 

and 
for 



tion (|n31), the formulae I^FT^ and (|5.23b|l for "Gl 



for Zf'^, 



the formula E 



Y^'^, the formula 



and the definitions |S!3Hll, and fOHjl . The 

result is 



• (U^ 



n 



M' 



E 



4(2fc + l)!!_^5 <f^> 



^ OO 

A ^BA 



i-il 
ki 



{2k + 21 + l)n Mt ^^SA<fKL 



k+l+2 



r 



BA 



k + 1 



'■Kr k+2 

"^BA 



4fc(2fc+l)!! ^g 



,BA 
<sK> 

r' 

' BA 



+ Oie^). (6.20) 



VII. CONCLUSION 

In this paper, we have given a surface integral deriva- 
tion of the full post-l-Newtonian DSX laws of motion 
(|4.3b|) . We have shown that these laws of motion ap- 
ply to a wide class of strongly self-gravitating objects, 
provided that the mass and current moments are appro- 
priately defined in terms of the asymptotic weak field 
metric in the buffer regions around each body. We have 
given an explicit form for the coupled equations of mo- 
tion of the bodies' center of mass worldlines including the 
effects of all the post-Newtonian mass and current mul- 
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tipole couplings. To the best of our knowledge this is the 
first time these equations of motion have been written 
out explicitly. The second paper in this series will in- 
clude a surface integral derivations of the evolution laws 
(|4.3all an d (I4.3c|l for the energy (mass monopole) and the 
spin Sf 
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APPENDIX A: USEFUL IDENTITIES 
INVOLVING STF TENSORS 

In this appendix we give the general definition of the 
STF projection of an arbitrary tensor. We also review 
some identities that are useful for manipulating expres- 
sions involving STF tensors. 

The STF projection of any tensor Tl is obtained by 
taking the symmetric part of T^,, and then subtracting 
out all the partial traces. One obtains in this way a 
unique symmetric tensor that is trace free on all pairs of 
indices. The general formula for this projection is 

[1/2] 
k=0 

(Al) 

where [1/2] is the largest integer less than or equal to 1/2, 
the coefficients are given by 

. ( 1)^ ^! (2/-2fc^l)!! 

^ ' {l-2k)l{2l-iy.l{2k)U' ^ ' 

III means 1{1-2){1^4) . . . (4)(2) or /(/-2)(/-4) . . . (3)(1), 
and 

Sl = T^L) (A3) 
is the symmetric part of T . For example, 

T^abO = Sabc — - [SabScdd + ^acSbdd + 5bcSadd\ ■ (A4) 
5 



From this peeling formula one can obtain the identities 

21 + i 



and 



21 + 1 



1 



-T,,i 



il + l){2l + l) 



iL 1 



(A6) 



(A7) 



which are valid for any STF tensor Ti+i. 

Next, some useful formulae involving derivatives are 



\x?dr 



1 



-{2l-l)d^L>\x\ 



and 



diL\x\ = d<^L>\x\ + ^[TpY^'^'^^^^R 



(A8) 



(A9) 



Also, given a sequence of STF tensors T^, one for each I, 
we have the identities 



1=0 



and 



1=0 '- ^ 

+TLd<jL>\x 



21 + 3 



(AlO) 



(All) 

These identities are used in Sec. IIII (Tl above in the com- 
putation of the transformation laws for the multipole and 
tidal moments. In Sec. IIVI we use the identities 



X X 



\l+2 



{aiT,L-in^-^ - biTLU^'^) , (A12) 



-.L-l 



(A13) 



where ai = {-1)H{21 - 1)!!, hi = (-1)'(2; + 1)!!, and 
n* = x'/la;!- We also use the following integrals over the 
unit sphere given in Thorne j38j | 



1 

47r 



(A14) 



From the definition of the STF projection one can de- 
rive the following "peeling formula" 2J for any STF ten- 
sor Tl and vector Vi: 



47r 



21 _|_ ^^{ni2---^i2l-li2l)^ 



(A15) 



21 



{l + l){2l + l)^''^'<'^-'^''^>^- 



- f T,S,n-n^ dn ^ (2/ + l)n 



II 



TlSl if fc = / 
if fc ^ (A16) 
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and 



= if |/c-Z| 7^ 1. (A17) 

Another useful peeling identity, valid for a A; + 1 index 
tensor which is STF on its last k indices, is 

T<:iKZL-K> — -J—^Ti^KZL-K> 
+ J-^T^K+iZL-{K+1)>Zi 



1l>l 



2{k + l){l-k) 
{l-k){l-k-l) 

(A18) 

Finally, for any STF tensor T and vectors Vi and we 
have the identity 



l + k + l 

fcq + i) ^ 

'{l + k+ 1){21 + 2k + i)^™^™^J^J<''i>^-2^ 



(A19) 



APPENDIX B: DERIVATION OF GAUGE 
TRANSFORMATION PARAMETERIZATION 

In this appendix we consider harmonic, conformally 
Cartesian coordinate systems on a spacetime region V x 
(to, ^i), where I? is a simply connected spatial region and 
{to,ti) is an open interval of time. We show that the 
most general gauge transformation between two such co- 
ordinate systems is of the form given by Eqs. (|2.17(l - 
(|2.18cll . up to constant displacements in time and up to 
time-independent spatial rotations. 

We start by reviewing the well-known argument that 
gives this result to Newtonian order. Let the coordinate 
transformation to zeroth order in e be 



t = t{t,x^) + 0{e^). 



(Bl) 



Substituting this into the metric expansion (|2.4|) , we find 
that the leading order expression for the spatial metric is 



1 dt dt 

dx^ dxi 



dx'dx' + 0(1). 



(B2) 



This is in conflict with the expansion H2.16|l unless 
dt/dx^ = 0. Similarly, the leading order expression for 
the time-time piece of the line element is 



dp 



0(1), 



(B3) 



which disagrees with the expansion H2.16|l unless dt/dt = 
±1. Assuming that the coordinate transformation pre- 
serves the time orientation and neglecting constant dis- 
placements in time we obtain t = t + 0{e'^). Therefore 
we can write 



t = t + e^a{t,x^) + 0{e'^), 



(B4) 



where the function a{t,x^) is as yet undetermined. 

The leading order expression for the spatial metric is 
now 



dx'' dx^ 
' dx"^ dxi 



dx'dx' + 0{e^) = Sijdx'dx^ + 0{e^), (B5) 



where we have used the expansion (|2.1t)|) . Thus, for each 
fixed the function a;* = (f, x^ ) is an isometry of 3- 
dimensional Euclidean space. It follows that 



= R- At)x^ + z\t) + 0(e2) 



(B6) 



for some time-dependent rotation matrix W ^it) and 
some time-dependent displacement z^(t). Using Eqs. 
H2.4(l , (|B4p and ljB6|l the leading order expression for the 
space-time piece of the line element is 



2<5,fci?^(t) iVAT)x^ +i\t) 



^^]didx'+0{e\ 
(B7) 

The first term here must vanish in order to be compatible 
with Eq. (|2.1()|) . which gives 



5,fci^^(^) tAt)x^ + zAT) 



da 

dx'- 



(B8) 



If i?' j (t) is non- vanishing, it is impossible to find any 
function a{t,x-') which satisfies this equation, since the 
left hand side is not a pure gradient. Therefore we con- 
clude that the rotation matrix is time-independent, and 
we choose the new coordinate system a;* so that i?' j — Sj . 
We can now solve Eq. IjBSp for the function a, which gives 



(B9) 



^{t,x^)^ac{t) + i\t)x\ 



where adt) is an arbitrary function of t, cf. Eq. (|2.18a|l 
above. 

To summarize, the coordinate transformation to New- 
tonian order is given by 

xHi,x^) = x' + zAt) + 0{e^) 
t{t,x^) = t + e^a{t,x^) + 0{s^), (BIO) 
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where a is given by Eq. (jB9|) . The transformation law for 
the Newtonian potential $ can now be obtained by sub- 
stituting Eqs. I)B10|I into the metric expansion (|2.4|l and 
comparing the time-time piece with the metric expansion 
(|2.16l) : the result is given by Eq. (|2.20a|l . 

We now turn to the post-Newtonian extension of this 
computation. We assume that the coordinate transfor- 
mation can be written as the Newtonian-order coordinate 
transformation IjBlOII plus arbitrary post-Newtonian cor- 
rection terms: 

x' = x' + z'{t) + e^h'{i,x^) + 0{e'^), 
t ^ t + e^a{t,x^)+e^f3{t,x^) + 0{e^). (Bll) 

Here the functions h^{t,x^) and /3{t,x^) are arbitrary. 
As before we can compute the transformed metric by 
combining the coordinate transformation (jBll|) with the 
metric expansion (|2.4|) . The resulting leading order ex- 
pression for the spatial metric is 



-2$(5,,; 



^dx^dx^, (B12) 



where we are using the notation H2.21|l . Comparing this 
with the metric expansion H2.16|l and using Eqs. H2.2Ua|) 
and H2.18a|) gives the following differential equation for 
h^(t,x^y. 



-25r 



ZkX 



— z 

2 



+ ZiZj. (B13) 



The general solution to this equation consists of a homo- 
geneous solution plus an inhomogeneous solution. The 
homogeneous solution is just the general Killing vector 
of three dimensional Euclidean space, which is 



h'{t, x^) ^hl{t) + e'^'^x^Rkit), 



(B14) 



where the functions hl(t) and Rk{t) are arbitrary. The 
full solution that we will use for /i* is the sum of (jB14|) 
and the inhomogeneous solution, which can be obtained 
by inspection. The result is [cf. Eq. 12.18b|l above] 

h\t,x^) = hiit) + e'^^XjRk{t) + ]^^{t)xjX^ 

-x'adt) - x'xjzdt) -I- ix'ij(<)i-'(<) 



+ -z\t)z^{t)x,. 



(B15) 



Next, we use Eqs. (^31, (ElSl), and l|5T5)l to 

compute the transformed gravitomagnetic potential CJ' . 
The result is 



J\ - 



1, 



it)ij{t)x' + -z'{t)XjX 



+x' [2zdt)zj{t) - x^Zj{t) - a^t)] 
~z\t) A^{t,x^) + 2adt} + ^x^Zjit) 



-P{t)z,it) 
+e^jkX^R''{t) 



§iit,x^)+m 

-e,jkR\t)z''{t). 



(B16) 



Combining this with the expression (|2.20a() for the trans- 
formed Newtonian potential, and using the harmonic 
gauge condition H2.6|l applied to both the original and 
barred coordinate systems gives the differential equation 



^^f3^-z,{t)x^+adt). 



(B17) 



The general solution to this equation is 



l^-zdt)x'' + ^adt) 



(B18) 

where /3h is an arbitrary harmonic function, cf. Eq. 
H2.18c|l above. This completes the derivation. 



APPENDIX C: PIECE OF SURFACE INTEGRAL 
THAT DEPEND LINEARLY ON MOMENTS 

In this appendix we compute explicitly the piece of 
the surface integral 14.4311 that depends linearly on the 
multipole, tidal and gauge moments. That linear piece 
appears on the right hand side of Eq. (|4.44|l as the func- 
tion Gii ''Ml, "Gl,Hl, Sl, 'Al, 'i^'l-.R)- We will show 
that this function vanishes. 

We start by noting that the splitting of the surface inte- 
gral H4.43|l into pieces that are linear in the moments and 
pieces that are quadratic in the moments is unambiguous 
for all the multipole and tidal moments, except for the 
Newtonian mass dipole "Mj(i). That mass dipole is con- 
strained by the Newtonian equation of motion (|4.2b|l . 
and therefore a term proportional to the fourth time 
derivative of "Mi could be re-expressed as a quadratic 
expression in the moments "Ml and "Gl and their time 
derivatives up to second order. We resolve this ambigu- 
ity by demanding that there be no dependence on "Mi 
in Qi ; the relevant term if present can be re-expressed as 
a quadratic expression and moved into the function J^i. 

To compute the linear piece of the surface integral, we 
simply drop all the quadratic source terms in the post-2- 
Newtonian field equations and gauge conditions H4.13|l - 
(HHni). We also drop the term p'T*-'' in Eq. lH^l^ . since 
the expression (|4.22|l for p^T*-' is explicitly quadratic. We 
also assume without loss of generality that "Mi = 0, for 
the reason discussed above. This yields from Eqs. 
and H4.44|l the set of equations 



1 

16^ 



9 3^0 



where 



djXo 



Co I ^iCo 



y'xo' = 0. 



(CI) 



(C2) 
(C3) 



Here the subscripts indicate that the post-l-Newtonian 
mass dipole associated with the potentials (<i>0 5 Coi V'o) 
vanishes, cf., the discussion in Sec. IIV El above. 



42 



To compute the function Gi^^ , we can pick any solution 
of the post-2-Newtonian equations HC2|I - HC3|I . since we 
showed in Sec. IIV El that the result is independent of 
which solution is chosen. A particular solution {^l,Xi) 
of Eqs. ljC3|) can be obtained using the expansion 13.5c() 
of the gravitomagnetic potential C'. This gives 



11 
Xi 



0, 

oo 

E 

1=0 



(C4) 



(-1) 



(+1 



2U 



2(2/ + 3);! 



(C5) 



These potentials do not satisfy the gauge conditions ljC2p . 
but we can fix this by adding appropriately chosen solu- 
tions of Laplace's equation. Thus, we define 



a = SI 



°° f-lV+i 1 1 

^-0 ■ II- 



Xo 



11 



E- 

i=0 



1 



1 



B,yLa;^(C7) 



Inserting Eqs. HC4p - (|C7|I into the formula IjCip gives 
1 1 1 • i?"^ • 

Qi — 4* — ^'t'jj — -^ijj- (^^) 

^ One mig ht tMik that tlti easiestl^ay to evaJt^te the expression 

ICll for Qi is to use Gauss' theorem to convert the surface inte- 
gral to a volume integral. However, this strategy does not work: 
Because the fields are only defined on the domain rg < r < ri one 
obtains a surface term at r = ro in addition to the volume term. 



To obtain the moments Xi, Cijj and Bijj we substitute 
Eqs. - l|('7p into the gauge conditions (|('2p . A useful 
intermediate result is 



E 

OQ 

E 

1=0 



2"MLdL\x\ 



2/ + 3^|a;| 



(2Z-|-3)/l 



[2\x\'^'^GL~i>iLX^] 



(C9) 



This gives Bijj = —Vi and Xi — Cijj/i — 0, yielding from 
Eq. ijUHl that Gi^Q. 



APPENDIX D: LAW OF MOTION FOR A 
SINGLE BODY WITH WEAK SELF-GRAVITY 



In this appendix we sketch briefly the DSX derivation 
O of Eq. translated into our notation. This 

equation is obtained by direct computation of the second 
time derivative of the mass dipole, making use of the 
stress-energy conservation law in the interior of the body. 
In our notation, the total mass dipole is [cf. Eqs. I|3.2|) 
and H3.24|l above] 

It is impossible to extend the definitions of the fields smoothly 
all the way to r = 0, so one is always forced to evaluate a surface 
term. 



J 



d?2 



6 df^ 5 dt 

The conservation equations V^T'^" can be written using the expansions (|2.4|l and (|2.2() in the form 



(Dl) 



at ^ ' ox^ at 



(D2) 



^ [(1 - 4e2$) ( + £2 P-T"')] + _^ [(1 _ 4£2$) ( _^ ^2 p^^j)] ^ 



These conservation equations can be used to evaluate explicitly the time derivatives appearing in Eq. HD1|I . Some 
algebra leads to the following expression 



d^'x 



^ > ax^ 2 ox^ 



(D4) 
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Taking two time derivatives of this expression and using the conservation equations HD2|I and ljD3|) gives 



UL J r<ir — L 



d 



d 



d 



d 



dx^ dx 



dt 



9$ 



4$T"' + x"T""— I (Fx 
at 



1 



(D5) 



We next substitute in exphcit expressions for the gravitational potentials in which the intrinsic terms are expressed 
in terms of integrals over the matter distribution using the field equations (|2.8a|l - H2.8c(l : 



$ = 
and 



\X — X' 



(D6) 



1=0 



2(2? + 3) 



°Gz 



LX 



(D7) 



It is a straightforward exercise to show that all the terms involving double integrals over x and x' in Eq. ljD5|l cancel 
out. The laws of motion can thus be obtained by simply substituting the tidal pieces of the gravitational potentials 
into Eq. (|D5|I . The remaining integrals over x can then be expressed in terms of the moments "Ml, ^"Ml and ZiL 
via the integral definitions H3.2|l . H3.24|l and 



Z,Lit)=4 'T°'{t,x')x<^>d^x. 



(D8) 



This gives 

oo ^ 



(=0 



(? + 2)(2Z + l)„ 



(2^ + 3) 



a + l')(2? + 3)'^^ + ^^.-^^[^1^ - Z^L '-Gl - "Gl 



M,L "GL-i2l + l) "Gl - I "M»L "Gl 



(D9) 



where we have used Eqs. (|4.2b|) . H4.46|l and (|D5|) . Using the STF decompositions (|3.15|l and (|3.1t)|) of the moments 



ZiL and YiL , it is straightforward to check that Eq. (|D9|) is equivalent to Eq. ()4.3b|) 



APPENDIX E: FORMULAE FOR MOMENTS IN 
TERMS OF SURFACE INTEGRALS 



s '^<i>^X' d'^^ = (^TTTW^*^ (2^ + 1)!!^^^' 



In this appendix we show that the various moments 

(E2) 

are uniquely defined by the expansions H3.5a|l - H3.5cl) , by ^^'^ 

writing down surface integrals from which the moments r 2 

can be explicitly computed. From the definitions H3.5a|l . f n^Lydjifi d = , . , ^ ^"Ml 

Kfi^ and (El we obtain ^ 

1 ■ IR'^' r ^ -1 



(; + l)(2Z + 3)i?''" (2^ + 1)!! 

^ 1 Ij^i+i l_l ■■ {I + 2)i^'+3 

9,$d2s, = + (2^ + 1)!!'^^' " 2(2/ - 1)(2; + 2(2/ + 3)!! 



(El) (E3) 
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Here is the natural surface element determined by 

the flat metric (dx^y + {dx^Y + {dx^Y, and the 2-surface 
S is the coordinate sphere r = R. By evaluating the right 
hand sides of these equations at several different values of 
R, and by using the decompositions l|3.15|l and (|3.16() . one 
can extract explicit expressions for the moments "Ml, 
"Gl [Eq. (Eil)], Hl, Sl, f^L, VL [Eq. dH], and p"Ml, 
^°Gl [Eq. IjESp ] in terms of the surface integrals and their 
time derivatives. 



If we substitute this equation into the full equation of 
motion H6.11(l , the following coefficients take the following 
new values 



APPENDIX F: COEFFICIENTS OF FINAL 
EQUATION OF MOTION AFTER 
SIMPLIFICATION USING SPIN EVOLUTION 
EQUATION 

In paper II j37|. the following spin evolution equation 
is derived 



B^A k=0 1=0 



n 



BA 

<mKL> 
„k+l+2 ■ 
' BA 



(Fl) 



i3)T)AB 



kW.M' 



(-l)''(2fc + 2^ + l )!! 



A kaB 



-2MfrM 



{21 + 3) 



4 ^B _c;aMI 



M 



( 4(1 - Sok) ^B 
fc + 2 '"^ 



M 



1 + 2 



^aMI 

m^A 



(2^2 + 3^ + 5) BA,-.A ,.B 



(F2) 



and 



~r,ABC _ (-l)'=+''(2fc + 2Z + l)!!(2p + 2g+l)! 



iKLPQ 



k\l\p\q\M' 



a + 2)(2/ + i) 

{21 + 3) 



A new three-body term is generated, which contributes to df^'^'. It can be written as 



(F3) 



oo ctt oo 



,BA 



\ \ \^ (7)'fsABC "■<iPQ> "'<jL> 
2^ 2^ 2^ ^JLPQ v+q+2 J+2 



1=0 p=0 q=0 



' CB 



' BA 



(F4) 



where 



(7)t,abc 



-l)P+i(2Z + l)!!(2p + 2g+l)! 



2MtMfQM%. 



(F5) 
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Finally, note that the 5^ terms in the second and last 
lines of Eq. (|6.13c|l cancel each other out. Therefore no 



three-body terms are generated from 
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TABLE II: In this table we list, for ease of reference, some of the symbols used in the paper in alphabetical order. We do not 
list symbols whose meaning is very conventional, or which are used only in the immediate vicinity of where they are introduced. 
For each item listed, we give a brief description, and also a reference to the equation or section in the text where the symbol 
first appears, or in the vicinity of which the symbol is first introduced. 



Symbol 



Meaning 



First appears in 



pn 
A 

a{t, X-') 
ac{t) 

Xij 

£ 

FUt) 

"GL{t) 

GUt) 
h\t,x') 

Kit) 
Mt) 

K 
L 

\L{t) 

A!(i) 

ML{t) 

MUt) 
N 

nf^(t) 

PL{t) 

p 

Ph 

$ 

V' 


Rk{t) 



Superscript appended to a symbol denoting that it is defined with respect to the global coordinate 
system 

Superscript prepended to a symbol denoting the Newtonian piece of a quantity 

Superscript prepended to a symbol denoting the post-l-Newtonian piece of a quantity 

Index appended to a symbol indicating that is associated with the A^^ body in an A'^-body system 

{B and C are used similarly) 

Function appearing in gauge transformation which parameterizes Newtonian-order changes in the 
time variable 

Piece of a{t,x-') that is independent of spatial coordinates 

Function appearing in gauge transformation which parameterizes post-l-Newtonian changes in the 
time variable 

Piece of /3{t,x-') that satisfies Laplace's equation 

symmetric tensor parameterizing the post-2-Newtonian spatial metric 

Post-Newtonian dimensionless expansion parameter 

Tidal moment of order I parameterizing tp about a worldline 

tensor density sometimes called the "gothic metric" equal to y/—gg^^ 

Newtonian gravitoelectric tidal moment of order I 

Post-l-Newtonian gravitoelectric tidal moment of order I 

Total gravitoelectric tidal moment of order I 

Free function in gauge transformation which parameterizes post-l-Newtonian translations 
Piece of h^{t,x-') that is independent of spatial coordinates 

Tensor density appearing in the Landau-Lifshitz formulation of general relativity 
Post-l-Newtonian gravitomagnetic tidal moment of order I 
Tidal moment of order I parameterizing tp about a worldUne 
The multi-index 61 62 ... 
The multi-index a\a2 ■ ■ .ai 

Intrinsic-type multipole moment parameterizing the harmonic gauge-transformation function 
Pi.{t,x^) 

Inertial moments that appear in the transformation law of '^Ghit), nonvanishing for I = 0,1 only 
Inertial moments that appear in the transformation law of YiL(t), nonvanishing for i = 1, 2, 3 only 
Inertial moments that appear in the transformation law of ^"GL{t), nonvanishing for i = 0, 1, 2 only 
Newtonian mass multipole moment of order I 
Post-l-Newtonian mass multipole moment of order I 
Total mass multipole moment of order I 

Mass multipole moment of order I defined in a body-frame non-rotating with respect to distant 
stars 

Intrinsic gauge moment of order / 
The multi-index aia2 . . . a„ 

Intrinsic multipole moment of order I parameterizing t/) about a worldline 

J**" component of a spatial unit vector pointing from the worldline of body A at time t to the 
worldline of body B at time t 
Tidal gauge moment of order I 
The multi-index ci C2 . . . Cp 
Momentum enclosed by a surface E 

Intrinsic multipole moment of order I parameterizing tp about a worldline 
Newtonian potential 

Post-l-Newtonian correction to the Newtonian potential 
The multi-index d\d2 . . .dq 

Function appearing in gauge transformation which parameterizes post-l-Newtonian rotations 



Sec.lTcl 



10 
(in 

Il2.18bll 

(EH 

Sec. EH 
Sec.lTFl 
llT?7ll 

lion 

Il3.36bll 

JO 

ll3^ 

Sec.im 
Sec.ira 

Sec.im 



Sec.im 
Il2.18bll 
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Symbol Meaning First appears in 

vba coordinate distance between the center of mass worldlines of bodies A and B, defined with respect Sec. II CI 
to the flat metric Sij 

sa time coordinate of a coordinate system adapted to body A 15.11 1 

Shit) Current multipole moment of order I 13. Set 

SL(t) Current multipole moment of order I deflned in a body- frame non-rotating with respect to distant 15.361 1 
stars 

t time coordinate of generic, harmonic, conformally Cartesian coordinate system in Sees. IITI- IIVI 12.41 

t time coordinate of global coordinate system for an A*' body system in Sees. IVI- IVll Sec. IV A1 

r*^" Components of the stress-energy tensor 12.211 

T'^" Landau-Lifshitz pseudotensor 14.41 

TL{t) Tidal- type multipole moment parameterizing the harmonic gauge-transformation function /3h(t, x-') 13.271 1 

T^{z) Taylor coefficients of the function \z — x\ about x = 15.181 1 

f// (t) Rotation matrix describing the dragging of asymptotic rest frames 11.31 1 

vf Velocity of the Ath body KTEi 

vf^ Relative velocity of bodies A and B 15.181 1 

a;* spatial coordinates of generic, harmonic, conformally Cartesian coordinate system in Sees. ITTl- liVI 12.41 

spatial coordinates of global coordinate system for an A'^ body system in Sees. IVl- IVll Sec. IV A1 

i^^ component of the post-2-Newtonian correction to the gravitomagnetic vector potential 14.111 1 

y-'^ spatial coordinates of a coordinate system adapted to body A Sec. II CI 

YiL (t) Tidal moments of order / of the gravitomagnetic potential 13.81 1 

Zi{t) Free function in gauge transformation which parameterizes Newtonian-order translations (12.171 1 

Zi{t) The center-of-mass worldline of a body, to Newtonian order 15.41 1 

'""zi{t) Center of mass worldline of a body, to post-l-Newtonian order Sec. IV d 

ZiL Intrinsic multipole moments of order I of the gravitomagnetic potential 13.81 1 

(^i i"^ component of the gravitomagnetic vector potential 12.41 1 



